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13  ABSTRACT 


Progress  is  reported  in  the  characterization  of  electronic  transport  and 
dielectric  breakdown  properties  of  technologically  important,  thin  insulating  films  on 
Si,  namely  SiOjj,  A$2®i*  SijN^  and  their  composites. 

New  corona-charging  techniques  have  been  developed  for  the  study  of  non-destruc¬ 
tive  breakdown  of  un-metallized  Si02  films,  enabling  the  measurement  of  I-V  character¬ 
istics  in  the  steady-state  discharge,  the  identification  of  the  current-carrier  in 
the  oxide,  an  accurate  study  of  electron  tunneling  from  the  Si  substrate,  and 
producing  strong  evidence  against  avalanching  in  the  oxide. 

Locally  destructive,  self-quenching  breakdown  of  AU-AI2O3-SI  and  Al-Si02~Si 
structures  has  baan  studied.  The  time  scale  for  breakdown  is  hundreds  of  nanoseconds 
for  the  Si  in  inversion,  and  tens  of  nanoseconds  in  accumulation.  Pre-breakdown  C-V 
(charge  storage)  measurements  suggest  breakdown  of  the  AI2O3  structure  at  the  Si 
interface,  and  of  the  Si02  structure  at  the  A1  interface.  A  silicon-oxide  interface 
under  a  thin  film  has  been  successfully  probed  in  the  SEM  with  high  lateral  resolution. 

Strong  electron-trapping  in  Al-implanted  Si02  films  has  been  verified  by  charge- 
discharge  techniques,  using  optical  charging  and  optical  or  thermal  discharging,  and 
has  been  successfully  modelled  theoretically.  C-V  measurements  have  confirmed  strong 
electron-trapping  in  Si02  films  previously  exposed  to  kilovolt  electron-beam  injection. 

Theoretical  modeling  of  self -quenching  breakdown  is  under  study  with  tunneling, 
joule-heating  and  their  interaction  (feedback)  as  the  key  ingredients.  The  dynamical 
behavior  of  runaway,  hot-electron  distributions  has  been  extensively  studied  by 
graphical  and  computer  techniques. 
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Abstract 


Progress  is  reported  in  the  characterization  of  electronic 
transport  and  dielectric  breakdown  properties  of  technologically 
important,  thin  insulating  films  on  Si,  namely  SiO_ ,  A1„0  Si  N, 
and  their  composites.  J  ^ 

New  corona-charging  techniques  have  been  developed  for 
the  study  of  non-destructive  breakdown  of  un-metallized  SiC^  films, 
enabling  the  measurement  of  I-V  characteristics  in  the  steady-state, 
the  identification  of  the  current-carrier  in  the  oxide,  and  the 
accurate  study  of  electron  tunneling  from  the  Si  substrate. 

Evidence  is  adduced  against  avalanching  in  the  oxide. 

Locally  desctructive ,  self-quenching  breakdown  of 
Au-A^O^-Si  and  Al-SiC>2-Si  structures  has  been  studied.  The 
time  scale  for  breakdown  is  hundreds  of  nanoseconds  for  the  Si  in 
inversion,  and  tens  of  nanoseconds  in  accumulation.  Pre-breakdown 
C-V  (charge  storage)  measurements  suggest  breakdown  of  the  AI2O3 
structure  at  the  Si  interface,  and  of  the  SiC>2  structure  at  the 
A1  interface.  A  silicon-oxide  interface  under  a  thin  film  has 
been  successfully  probed  in  the  SEM  with  high  lateral  resolution. 

Strong  electron  trapping  in  Al-implanted  Si02  films 
has  been  verified  by  charge-discharge  techniques,  using  optical 
injection  of  charge  carriers  into  the  films  for  charging  of  the 
traps  and  optical  or  thermal  discharge  of  the  trapped  carriers. 

The  trapping  has  been  successfully  modeled  theoretically. 

Theoretical  modeling  of  self-quenching  breakdown  is 
under  study  with  tunneling,  joule-heating  and  their  interaction 
(feedback)  as  the  key  ingredients.  The  dynamical  behavior  of 
runaway,  hot-electron  distributions  has  been  extensively  studied 
by  graphical  and  computer  techniques. 
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I .  Int reduction 

This  report  details  progress  in  a  comprehensive 
experimental  and  theoretical  research  program  for  the  character¬ 
ization  of  electronic  transport  and  dielectric  breakdown 
properties  of  thin  insulating  films  on  silicon  which  are 
technologically  important  (SiO^,  Al^O  ,  Si^N^,  and  their 
composites ) . 

The  experimental  approaches  being  used  are: 

i)  Completely  nen-des4- ructive  breakdown  of  the  insulating  film 
achieved  by  contacting  the  exposed  surface  of  the  film  by  a 
corona  discharge  in  a  suitable  gas, 

ii)  locally  destructive,  self-quenching  breakdown  of  the  film, 
achieved  with  thin  metallization, 

lii  lopticai  injection  of  charge  into  the  film  and  optical  or 
thermal  discharge  of  trapped  carriers  in  the  film, 

iv )  electron-beam  injection  of  kilovolt  electrons  stopped  inside 
the  film  under  a  thin  metallization, 

v)  elemental  analysis  and  topographic  imaging  of  breakdown- 
induced  damage  in  the  film, 

vi  1  electron-beam  studies  of  the  solid-solid  interface  between 
the  silicon  substrate  and  the  film, 

The  theoretical  studies  which  are  being  carried  out  are: 

1)  Monte-Carlo  calculations  of  hot-electron  distributions 
produced  by  high  electric  fields  in  the  film, 
ii )  modeling  studies  directed  toward  identification  of  those 
physical  mechanisms  critical  to  the  electrical  breakdown  of  the  film, 
lii ) tunneling  of  electrons  out  of  a  triangular  potential  well 
in  the  Si,  produced  by  strong  band-bending  at  the  interface  into 
the  film. 

Highlights  of  the  current  report  are: 

A.  Mon-Destructive  Breakdown 

In  the  non-destructive  breakdown  of  SiQ^  films  thermally 
grown  on  Si,  the  sign  of  the  current-carrier  in  the  oxide  has 
been  determined  by  use  of  a  shallow  P-M  junction  diode  as  the 
substrate.  For  both  positive  and  negative  coronas,  electrons 
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s re  the  current  carrier.  With  negative  corona  it  has  been 
directly  verified  that  each  electron  entering  the  silicon  from 
the  cxide  genera+es  approximately  one  electron-hole  pair  in 
thermalizing  to  the  bottom  of  the  silicon  conduction  band. 

A  direct-comparison  technique  has  been  developed  for 
measuring  the  surface  potential  of  the  exposed  film  under  steady— 
sfa+e  charging  by  the  corona-  With  this  technique  it  has  been 
shown  that  there  are  no  significant,  differences  in  the  I-V 
'-haracteristics  for  n-and  p-type  substrates. 

The  vibrating  Kelvin-probe  technique  for  measuring 
surface-potential  decay  rates,  after  charging,  has  been  developed 
mtc  an  accurate  and  reliable  vehicle  for  studying  tunneling  of 
ele  .crons  from  Si  into  SiOg.  This  technique  yields  sin  electron 
effective  mass  in  the  oxide  of  m^/in  =  0.39,  in  agreement  with 
previously  reported  values. 

Under  negative  corona  charging,  at  fields  above 
l.L  x  10  V/ cm,  positive  charge  is  found  in  the  oxide  through 
G-V  measurements  with  a  mercury  probe.  Arguments  are  adduced 
that  the  positive  charge  is  the  result  of  injection  of  hot 
holes  from  the  Si  into  the  SiOg,  the  holes  being  created  by  the 
incoming  electron  (the  primary  current-carrier  in  the  oxide). 

Two  independent  experiments  negate  avalanching  in  the 
Sx02  film  as  a  likely  significant  factor  in  local  film  breakdown. 
The  shallow- j unction  diode  experiment  with  negative  corona 
es+ablishes  that  a  significant  fraction  of  electrons  in  the 
oxide  cannot  have  as  much  as  even  2  ev  excess  energy.  Further, 
the  complete  absence  of  oxide  charging  with  positive  corona 
establishes  that  a  negligible  number  of  electron-hole  pairs 
can  have  been  created  in  the  breakdown. 

B*  Locally  Destructive,  Self-Quenching  Breakdown 

Locally  destructive,  self-quenching  breakdown  has  now 
been  studied  on  Au-Al^O^-Si  and  Al-SiO^-Si  structures,  where  the 
A-*-2C3  py  ”0  lytic  ally  deposited  and  the  SiOg  thermally  grown. 

A  result  common  to  both  structures  is  an  order-of-magnitude 
difference  in  the  time  scales  for  breakdown  between  those 


situations  where  the  Si  substrate  is  in  accumulation  at  the 
interface  (tens  of  nanoseconds)  and  where  the  substrate  is  in 
inversion  (hundreds  of  nanoseconds).  The  corresponding 
equivalent  resistances  in  the  discharge  path  are  roughly  one 
hundred  ohms  (Si  in  accumulation)  and  a  few  kilohms  (Si  in 
inversion).  In  the  case  of  AlgO^j  at  high  fields  at  pre¬ 
breakdown,  positive  gate  (Au)  voltage  charges  the  film  negatively 
and  negative  gate  voltage  charges  the  film  positively ,  as 
determined  by  C-V  measurements.  These  results  are  consistent 
with  breakdown  always  initiated  at  the  silicon  interface.  The 
results  for  SiO^  are  the  opposite:  at  high  fields  at  pre¬ 
breakdown,  positive  gate  (Al)  voltage  charges  the  film  positively 
(up  to  lCr  /cin  volume  density  of  charge)  and  negative  gate 
voltage  produces  relatively  little  charging  (C-V  flat-band  shift 
less  than  5  volts).  These  results  are  consistent  with  breakdown 
always  initiated  at  the  Al  interface.  Electrons  are  not  strongly 
trapped  in  a  clean  SiO^  film  and  therefore  do  not  charge  up  the 
film,  whereas  holes  are  always  strongly  trapped  at  the  Si-SiO^ 

interface  and  therefore  charge  up  the  film.  With  positive  gate 
bias  the  breakdown  was  sensitive  to  scratches  on  the  gate;  with 
negative  gate  bias  it  was  insensitive.  With  positive  gate 
bias,  optical  and  SEM  photography  of  the  exposed  breakdown  area 
revealed  it  to  be  round,  smooth-edged,  with  %  20u  diameter, 
and  a  distinct  inner  region  0.%  diameter;  further  the 
overall  size  of  the  breakdown  region  appears  to  be  independent 
of  whether  the  silicon  substrate  is  n-typeor  p-type .  With 
negative  gate  bias  the  exposed  breakdown  area  appears  round, 
rough-edged,  with  ^  5-<a.  diameter,  and  a  distinct  inner  region 
■u  l.2u  diameter.  Only  p-type  substrates  have  to-date  been 
studied  in  this  polarity.  The  exposed  area  of  breakdown  has 
been  probed  for  elemental  composition  in  the  SEM  using  x-ray 
fluorescence . 

C.  Secondary-Electron  Spectroscopy 

o 

With  relatively  thin  films  (<  1000A  thickness )  the 


SEM  has  been  used  to  image  the  silicon-oxide- interface .  The 
Te:hnique  employed  sweeps  the  detection  volume  through  the 
interface  region,  comparing  the  absorbed  current  images 
during  the  sweep.  Up  to  the  present,  lateral  resolution 
of  200A  has  been  achieved,  and  depth  resolution  sufficient 
to  image  detail  at  the  interface  but  remaining  to  be  determined 
quantitatively , 

D.  Charge-Discharge  Studies 

Electron  traps  in  aluminum-implanted- SiO^  films  have 
been  studied  by  electron  injection  induced  by  internal  photo¬ 
emission,  followed  by  optical  or  thermal  de-trapping.  The 
optical  de-trapping  required  hv  >  k.6ev.  Thermal  de-trapping 
was  achieved  at  temperatures  below  350°C.  Thermal  annealing  at 
600°C  substantially  reduced  the  concentration  of  electron  traps, 
suggesting  that  the  traps  are  located  at  defect  centers  created 
by  the  implantation  process.  A  model  of  electron  trapping  has 
been  formulated  which  illustrates  the  dynamic  features  of  the 
trapping  process  and  correctly  predicts  the  time  dependence 
of  the  flat-band  voltage  and  photocurrent  during  electron 
photo-injection;  both  analytical  and  numerical  calculations  have 
been  made  Application  of  the  model  to  the  experimental  results 
indicates  that  all  photomjected  electrons  are  trapped  in  the 
oxide  Also,  the  analysis  yields  a  numerical  value  for  the 
centroid  of  the  resulting  space-charge  distribution. 

E  Electron-Beam  Injection  Studies 

A  square-law  portion  of  the  steady-state  I-V 
characteristic  obtained  in  the  electron-beam  injection  experiment 
has  previously  been  interpreted  in  terms  of  electron-trapping 
in  the  oxide  in  equilibrium  with  photon  de-trapping,  the  photons 
being  generated  by  the  incident  (primary)  electron  beam  in  stopping. 
The  presence  of  electron  traps  in  the  electron-bombarded  oxide 
has  now  been  independently  confirmed  by  internal  photoemission 
experiments.  C-V  measurements  have  yielded  flat-band  voltage 
shifts  on  the  order  of  25  volts  corresponding  to  negative  charge 
in  the  oxide.  The  electron  traps  are  not  present  in  oxide  films 


which  have  not  been  bombarded,  by  the  kilovolt  electron  beam. 
Further  work  will  concentrate  on  obtaining  a  quantitative  picture 
of  the  production  of  traps  by  the  beam.  Toward  this  end,  an 
improved,  more  uniform  beam,  available  from  a  tungsten  cathode, 
is  under  design, 

F .  Theoretical  Modeling  of  Local  Breakdown 

On  the  theoretical  side,  the  major  effort  has  been 
spent  on  identifying  those  physical  mechanisms  which  appear  to 
be  basic  for  electrical  breakdown  of  the  films.  Two  ingredients 
are  essential  for  any  model  of  breakdown:  l)  injection  of 
carriers  into  the  film  from  the  substrate  or  gate,  and  2)  joule 
heating  in  the  film.  Generally,  models  of  breakdown  attempt 
also  to  identify  positive  feedback  mechanisms.  In  SiO^, 
breakdown  is  certainly  initiated  by  Fowler-Nordheim  tunneling 
into  the  film  from  the  Substrate  or  gate,  though  at  a  later 
stage  thermionic  emission  may  be  significant.  It  is  likely 
that  the  significant  tunneling  occurs  at  sites  of  irregular 
structure  (deviation  from  flatness)  at  the  cathode  surface 
where  there  is  field  enhancement,  A  simple  statistical  model 
is  given  for  estimating  the  distribution  of  protuberance  sizes. 
Current  initially  flows  in  filaments  in  the  film  inside  which 
there  is  rapid  heating  to  exceedingly  high  temperatures. 

(Because  of  the  small  filament  radius,  radiative  transfer  of 
heat  out  of  the  filament  is  insignificant  and  the  filament 
temperature  can  get  very  high.)  The  spreading  resistance  of  the 
silicon  substrate  is  an  important  element  in  the  evolution  of  the 
breakdown  in  that,  at  some  current  level,  it  brakes,  the  runaway 

7 

in  the  filament.  The  fields  are  high  enough  (:^  1  x  10  V/cm) 
that  ionic  conductivity  (e.g.,  of  Na  ions)  may  be  appreciable 
even  at  room  temperature.  Feedback  mechanisms  which  have  been 
considered  are:  l)  thermal  runaway  in  the  silicon  substrate, 
which  may  entirely  reverse  the  braking  action  of  the  spreading 
resistance  of  the  silicon;  2)  a  thermionic  instability  in  which 
the  thermionic  current  from  the  cathode  is  high  enough  to 
sustain  the  cathode  at  high  temperature,  and  which  provides  a 
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mechanism  for  substantial  radial  growth  of  a  filament;  3)  growth 
of  a  conducting  protuberance  with  consequent  increase  of  the 
field  enhancement  which  was  initially  determined,  at  room 
temperature,  by  the  pure  topography;  4 '  ion-induced  field 
enhancement  following  ion  motion  in  the  high  fields  at  elevated 
temperatures-  This  effect  has  recently  been  demonstrated 
experimentally.  It  is  a  very  powerful  mechanism  in  lowering  the 
breakdown  field  Impact  ionization,  historically  a  very  popular 
mechanism  in  the  theoretical  literature,  is  shown  to  be  highly 
unlikely  in  large  bandgap  insulators.  This  mechanism  is  also 
net  supported  by  the  available  experiments  (see  A). 

G.  Tunneling  Calculations 

Theoretical  work  has  aiso  been  done,  in  support  of 
th<=  non-destructive,  corona-induced  breakdown  experiments,  on  the 
prcblem  of  tunneling  of  electrons  out  of  the  triangular  potential 
well  produced  inside  the  silicon  at  the  Si-SiO^  interface  by 
strong  band-bending.  This  problem,  which  is  an  extension  of  the 
previously  solved  problem  of  space  quantization  inside  the 
triangular  well,  also  has  a  complete  analytical  solution,  in 
terms  cf  Airy  functions-  The  calculations  indicate  that 
tunneling  is  mostly  out  of  the  lowest  sub-band  of  energy  levels  in 
the  well. 

H.  Hot-Electron ,  Monte-Carlo  Calculations 

Seme  further  hot-electron,  Monte-Carlo  calculations 
have  been  carried  out  with  available  programs.  The  scaling 
law  for  the  constant  mean-free-path.  problem  has  been  more 
extensively  studied  for  large  scattering  anisotropies  in  the 
cne-dimensionai  random-walk  simulation.  Also,  a  more  systematic 
exploration  of  the  distribution  runaway  associated  with  energy- 
dependent  mean-free-paths  has  been  carried  out.  (Only  rudimentary 
results  on  this  problem  were  available  in  the  previous  report.) 

In  terms  of  a  simple  graphical  analysis  we  now  have  a  good 
picture  of  how  the  distribution  function  runs  away. 
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2 .  Non-Destructive  Breakdown  and  Transport  in  Un-Met  alii  zed  Films 
(Z.  Weinberg  collaborating) 

A.  Background  Discussion  and  Summary 

Recent  progress  has  been  successful  in  firmly 
establishing  the  corona  discharge  as  a  powerful  research  tool 
for  the  study  of  thin  insulating  films «  This  discharge  provides 
a  simple  means  for  charging  ,  positively  or  negatively,  the 
exposed,  unmetallized  surface  of  the  film.  The  film  is  easily 
charged  up  to  breakdown,  which  is  totally  non-destructive 
because  of  the  absence  of  lateral  conductivity  along  the  film 
surface.  Thus,  true  steady-state  and  quasi-steady-state 
currents  can  be  drawn  through  the  film  (the  latter  if  there  are 
slow  charging  effects  in  the  film).  There  are  several  unique 

i 

advantages  to  the  corona,  discharge  technique:  1)  the  experiments 
are  done  at  approximately  atmospheric  pressure  so  that  the 
charging  of  the  surface  takes  place  at  negligible  ion  kinetic 
energies  (  <  1  ev);  thus  damage  to  the  film  produced  by  ion 
kinetic  energy"'"  is  avoided,  2)  film  breakdown  is  reached  at 

1  x  10  V/cm,  corresponding  roughly  to  a  surface  ion  density 
13  2 

of  10  ,  cm  .  Thus,  the  effective  gate  associated  with  the  ion 
deposition  is  at  most  a  few  percent  of  one  monolayer,  and 
3)  because  of  the  absence  of  lateral  conductivity  along  the  film 
surface,  'weak'  spots  play  a  negligible  role  in  the  (non¬ 
destructive)  breakdown  of  the  film.  If  the  film  is  metallized, 
the  stored  energy  in  the  film  capacitor  discharges  through  a 
'weak'  spot,  producing  severe  local  destruction. 

The  results  obtained  to-date  have  established  that 
the  non-destructive  breakdown  of  SiO^  films  induced  by  either 
positive  or  negative  coronas  in  dry  air  is  due  to  the  tunneling 
of  electrons  from  either  the  silicon  (positive  corona)  or  from 
surface  electron  traps  (negative  corona)  into  the  film.  This 
has  been  established  by  use  of  a  shallow-,) unction  P-N  diode  for 
the  substrate.  The  breakdown  field  strength  of  the  film  is 
constant  over  the  range  of  thicknesses  studied  to-date,  namely- 
600  to  2700A,  and  the  value  of  the  breakdown  field  is 
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6  x  10  V/ cm.  for  positive  corona  (electrons  tunneling  from  the 

7 

silicon)  and  1.3  x  10  V/cm.  for  negative  corona  (electrons 
tunneling  cut  of  surface  traps).  With  the  positive  corona  an 
unusually  accurate  study  of  tunneling  from  the  silicon  into  the 
oxide  has  been  possible  using  a  vibrating  Kelvin  probe  to  monitor 
the  decaying  surface  potential  (after  corona  charging).  With 
negati  e  corona,  C-V  measurements  with  a  mercury  probe  have 
established  the  build-up  of  positive  charge  in  the  oxide  near 
the  interface.  At  present,  the  available  data  on  this  point  to 
the  tunneling  of  hot  holes  from  the  silicon,  these  holes  being 
created  by  electron-hoxe  pair  generation  by  the  incoming 
electrons  during  their  thermalization  (these  electrons  enter 
the  silicon  from  the  oxide  with  at  least  3  ev  kinetic  energy), 

B.  The  Corona  Discharge 

A  corona  discharge  in  a  gas  ambient  is  produced  by 

field  intensification  at  sharp  points  or  thin  wires  under 

application  of  high  voltage.  It  is  a  well-known  phenomenon  in 

the  electric  power  industry,  in  that  the  corona  generated 

around  a  high-voltage  transmission  line  is  an  undesirable  loss 

mechanism  -  in  fact,  the  corona  losses  determine  the  upper  limit 

cf  voltage  usable  for  electric  power  transmission.  More  recently 

the  corona  discharge  has  played  a  more  benevolent  role  in 

technology;  it  is  a  standard  component  in  electro-copying 

equipment,  being  used  to  charge  up  the  photo-sensitive  paper. 

In  general,  the  corona  is  a  convenient  source  of  electrostatic 

charge  in  a  gas  ambient.  It  has  been,  and  is  being,  used  by 
2 

R.  Williams  for  the  study  of  SiO^  films,  but  not  as  far-reachingly 

as  in  the  present  work. 

3 

Shahin  has  carried  out  a  mass -spectroscopic  analysis 
of  the  ionic  constituents  produced  by  a  corona  discharge  in 
room  air  at  various  humidities.  He  found  that  the  dominant 
ions  are  the  following: 

For  positive  corona, 

(H20)nH+  with  n  =  k  ,  5  ,  6 

Estimated  ionization  energy:  3.8,  3.2  ,  2.6  ev 
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For  negative  corona,  CO^  ,  with  en  estimated  electron  affinity 
of  3.8  +  1  ev. 

The  above  positive  ions  are  far  too  shallow  to  inject 
holes  into  the  valence  band  of  SiO^,  whose  maximum  lies  approx¬ 
imately  10  ev  below  the  vacuum.  {Only  a  handful  of  positive 

ions  are  deep  enough  to  inject  holes  directly  into  the  valence 

+  +  +  +  x  , 

band  of  SiO^;  these  include  He  ,Ar  ,Ne  ,Ng  ,  etc.;.  Likewise, 
the  above  negative  ions  are  far  too  deep  to  inject  electrons 
into  the  conduction  band  of  SiO^,  whose  minimum  lies  approximately 
1  ev  below  the  vacuum.  Thus,  ions  of  either  sign  can  only  charge 
up  the  surface  of  the  exposed  oxide,  until  non-destructive 
breakdown  sets  in.  In  the  terminology  of  energy  bands,  these  ions 
constitute  surface  states  lying  inside  the  oxide  bandgap.  Very 
high  fields,  ?  5  x  10^  V/cm.  are  required  to  allow  charge  to 
tunnel  from  these  states  into  the  corresponding  band  of  the 
oxide . 

C.  The  Basic  Electronic  Processes 

The  basic  electronic  processes  by  which  conduction 
current  can  flow  through  the  oxide  in  the  corona  experiments 
are  schematically  illustrated  in  Fig.  1(a)  and  (b).  With  the 
corona  positive,  Fig.  1(a),  the  silicon  substrate  will  be  in 
electron  degeneracy  at  the  interface  because  of  the  very  high 
field.  Thus,  n-type  Si  will  be  in  heavy  accumulation  at  the 
interface,  and  p-type  Si  will  be  in  deep  inversion.  Possible 
injection  modes  are:  l)  electrons  tunneling  into  the  oxide 
conduction  band  from  the  degenerate  electron  gas  in  the  conduction 
band  of  the  silicon  at  the  interface,  2)  electrons  tunneling 
into  the  oxide  conduction  band  from  the  valence  band  of  the 
silicon,  or  3)  holes  tunneling  from  the  positively  charged 
surface  states  (adsorbed  positive  ions)  at  the  exposed  surface 
into  the  valence  band  of  the  oxide.  At  any  one  field,  one  of 
these  three  injection  modes  will  be  dominant.  At  a  higher 
field  another  mode  may  take  over.  (Note  that  for  the  sake  of 
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simplicity  we  have  drawn  a  single  energy. .level,  of  positively- 
charged  surface  states;  in  practice  there  may  be  several 
energy  levels  corresponding  to  different  adsorbed  ions,  and 
further,  each  '  single'  level  may  in  fact  be  somewhat 
smeared  out  into  a  surface  band  due  to  various  interaction 
effects.  At  present  very  little  is  known  about  adsorbed  ionic 
levels  cn  the  surface  of  SiO^,  Figure  1(b),  corresponding  to 
negative  corona,  is  the  dual  of  Pig.  1(a).  The  silicon 
substrate  will  now  be  in  hole  degeneracy  at  the  interface  at 
the  high  fields.  Hence,  p-type  Si  will  be  in  heavy  accumulation 
and  n-type  Si  in  deep  inversion  at  the  interface .  Holes  may- 
tunnel  from  either  the  valence  band  or  the  conduction  band  of 
the  silicon,  or  electrons  may  tunnel  from  the  negatively 
charged  surface  states  (adsorbed  negative  ions). 

From  Figs.  1(a)  and  (b).it  is  clear  that,  with. either 
polarity  cf  corona,  either  sign  of .current-carrier  may  be 
dominant.  A  major  result  of  the  present  work  is,  in  fact,  the 
determination  of  the  sign  of  current-carrier  -.negative,  for  both 
polarities,  foi  corona  in  dry  air. . 

D.  The  Basic  Experiment 

A  schematic  of  the  basic  corona  experiment  is  shown 
in  Fig.  1(c).  The  ions  are  produced  in  the  immediate  vici.iity 
of  the  corona  needle.  The  dominant  element  in  the . subsequent . 
fate  of  the  ions  is  their  short-  mean  free  path  ).  for  collisions 
with  the  neutral  gas  molecules ;  namely  a  is  on  the  order  of 
10  ^  cm.  at  atmospheric  pressure.  As  a  result  the  ions,  remain 
thermalized  at  the  ambient  temperature,  corresponding  to.  a.. kinetic 
energy  =  0.025  ev.  This  very  low  kinetic  energy  of.  impacting, 
ions  protects  the  oxide  film  from  undesirable  excitations  which 
plague  other  methods  of  depositing -charge. on  the  film.  A  further 
consequence  of  the  strongly  collision^dominated.  motion  is  that 
the  ions  closely  follow  the  electric  field  lines.  This 
feature  has  provided  the  basis  for  measurement  of  the  surface  . 
potential  of  the  film  during  steady-state  charging  by  the 
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ccrena,  namely  by  a  comparison  method  described,  below. 

The  corona-induced  flux  of  ions -serves  as  an.  ideal 
current -limited  source.  The  exposed  film  surface  serves,  as  a 
perfect  floating  electrode  which  self-adjusts  its  potential, 
in  the  steady 'State, to  accept  the  ion  current  which  the  bulk 
of  the  film  can  transport  to  the  silicon  substrate. 

In  a  dry  ambient  there  is  negligible,  lateral 
conductivity  along  the  exposed  surface  of  the  oxide.  Thus, 

'weak'  spots  in  the  film  are  unable. to  discharge,  the.  electrostatic 
energy  stored  throughout  the  film.  Their  sole  effect  is  to 
create  a  highly  localized,  very  small  area  of  potential  deviation 
from  the  average  potential  over  the  film  surface.  The  total 
area  of  potential  deviation  is  obviously  proportional  to  the 
area  cf  the  weak  spots;  if  the  latter  is  small  comparted  to  the- 
total  exposed  film  area,  the  overall  effect  of  the  weak  spots 
is  correspondingly  small.  This  unusual  feature  of  the  corona 
method  is  in  startling  contrast  to  the  conventional  MQS. 
capacitor  studies,  where  a  single  weak. spot  will  lead  to. 
locally  destructive  breakdown,  because  of  the  high. lateral, 
conductivity  of  even  the  thinnest  metallization. 

The  samples  studied  to-date  have  all  been  SiO^  films 
thermally  grown,  in  a  dry  ambient,  on  (100)  waiters  of  silicon, 
both  n-  and  p-type  in  the  resistivity  range  1-2  ftem. 

Film  thicknesses  studied  were 

o 

p-substrate  :  501,  1027,  2570A 

n-substrate  :  587,  1000,-  2580A  ,- 

as  measured  by  ellipsometry  to  an  accuracy  better  than  1% .  The 
index  of  refraction  of  the  oxides  is  1.46.  Contact  to  the  back 
of  each  wafer  was  made  by  A1  evaporation  following-  the  appropriate 
shallow-dopant  diffusion  (boron  into  p-Si,  phosphorus  into  n-Si) 
to  provide  a  high-conductivity  region  of. contact.  The  sample 
dimensions  were  7.5  x  7.5  mm.  All  of  the.  samples,  were  provided 
by  the  RCA  laboratories  in  Princeton. 

All  experiments  were  done  in  an  environmental  chamber 
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which  allows  control  of  the  gaseous  ambient.  The.  chamber  is 
equipped  with  an  air  lock  for  transferring-material,  in.  and  out, 
and  a  pair  of  gloves  for  exterior  manipulation  of  the  interior 
equipment.  The  humidity  inside  the  chamber  is  kept  low  by 
circulating  the  interior  ambient  through  a.  desiccating  system 
containing  both  silica-gel  and  a  molecular  sieve.  The 
equipment  includes  provision  for  introducing  a  slow,  controlled 
flow  of  dry  Ng  or  ether  gases.  The  unit  can  be  completely 
covered  to  shield  the  experiment  from  exposure  to  ambient  light, 
if  necessary. 

A  schematic  of  the  corona  source  is  shown  in  Fi g ,  2. - 
The  source  of  the  discharge  is  a  gold  plated  needle  of. the  type 
used  as  a  microprobe  for  integrated  circuitry.  It  proved  to 
be  a  reliable  .and  stable  source,  and  provided  better  uniformity 
of  charging  ever  the  small  sample  area  than  a  thin  wire. 

The  surface  potential  of  the  exposed  film  was  measured 
under  two  conditions:  transient  decay,  following  charging,  and 
steady-state  charging.  The  transient  decay  was .  monitex-ed  by  a 
vibrating  Kelvin  probe,  as  illustrated  in  Fig.  3(a).  The 
vibrating  metal  plate  generates  an  ac  signal  proportional  to  the 
difference  in  potential  between  the  plate  (known)  and  the  film 
surface  (unknown).  A  lock-in  detector,  sensitive  only  to  the 
400  hz  frequency,  monitors  the  ac  signal.  A  precision  voltage 
supply  varies  the  potential  of  the  vibrating  plate  until  a  null 
is  obtained  in  the  ac  detector.  The  detection  sensitivity,  is 
+  0.05  V  relative  to  an  average  film-surface  potential  of 
100  volts. 

E .  Steady-State  Potential  Measurements 

The  steady-state  charging  potential-  is  measured  by  a 
comparison  technique,  as  illustrated  in  Fig.  3(b).  .The  sample 
and  a  metal  reference  plate  of  equal  area  are  placed  symmetric  ally 
in  the  corona  field  and  the  currents  to  both  are  monitored.  If 
the  surface  potentials  of  the  film  and  the  metal  plate  are  unequal 
then  the  two  currents  will 1  also  be-  Unequal-/  Conversely-/  -l: '■ 

equal  currents  correspond  to  equal  potentials.  In  the  imple- 
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mentation  of  this  simple  idea,  the  two  surfaces  are  surrounded 
by  a  grounded  mask,  as  shown  in  Fig.  3(b).  The  mask  ensures 
equal  exposed  areas  of  the  two  surfaces  and  prevents  ions  from 
reaching  the  sides  or  back  of  the  sample.  The  metal  reference 
plate  is  also  grounded,  the  current  being  controlled  by  a  grid, 
as  shown  The  potential  V_.  of  the  silicon  back  plate  is 

D 

adjusted  so  that  the  sample  current  Ig  is  the  sene  as  the  current 
I  collected  by  the  metal  plate.  At  this  point  the  exposed 
surface  of  the  film  must  also  be  at  ground  potential,  and  the 
oxide  voltage  V  must  cancel  the  applied  voltage  V :  V  =  -  V  . 

OX  Jj  OX  D 

In  practice,  the  sensitivity  of  this  technique  can  be  as  good 
as  +  IV.  However,  when  the  sample  response  is  very  slow,  the 
uncertainty  in  the  final  balance  hinders  its  reliability.  In  a 
regime  of  quick  response  the  method  is  both  accurate  and  con¬ 
venient.  In  this  respect  the  method  is  quite  complementary 
to  the  vibrating  probe  method,  which  is  rendered  impractical 
by  quick  decay.  In  fact,  development  of  the  steady-state, 
comparative  method  was  spurred  by  the  problem  of  coping  with 
deep-depletion  in  the  silicon  substrate.  With  the  substrate  in 
deep  inversion,  a  sufficiently  fast  rate  of  tunneling  of 
carriers  out  of  the  inversion  layer  into  the  oxide  can  drive 
the  substrate  into  a  steady-state  deep  depletion.  Since  a 
sizable  voltage  can  be  absorbed  across  the  deep-depletion  layer 
the  measured  total  sample  voltage  and  the  oxide  voltage  can 
differ  substantially,  leading  to  large  uncertainties  in  the  latter 
quantity.  Since  deep-depletion  in  silicon  decays  to  normal 
inversion  on  a  time  scale  shorter  than  a  second  at  room 
temperature,  the  Kelvin  probe  method  would  be  ineffective  in 
studying  this  important  regime.  The  steady-state,  comparative 
method,  in  contrast,  is  well  suited  to  this  purpose. 

F.  C-V  Measurements 

A  standard  technique  for  measuring  the  sign  and 
quantity  of  charge  stored  in  an  insulating  film  is  the  C-V 
measurement,  conveniently  made  at  a  frequency  of  1  Mhz. 


Contact  tc  the  insulator  surface  is  made  by  a  probe  in  the  form 
of  a  small  mercury  ball  held  in  a  platinum  container  with  a 
20  mil  diameter  opening.  The  physical  mechanics  of  the 
measurement  are  such  that  it  can  only  be  carried  out  after 
charging,  i.e.  during  the  decay  period.  As  is  well-knuwn,  the 
C-V  measurement  is  most  sensitive  to  charge  stored  at  or  near 
the  substrate-film  interface,  and  is  least  sensitive  to  charge 
stored  at  or  near  the  exposed  film  surface .  Note  that  the  sign 
of  the  charge  stored  in  the  film  after  corona  charging  need  not 
be  the  same  as  the  sign  of  the  dominant  current-carrier  either 
during  charging  or  decay.  Thus,  it  is  a  matter  of  common 
experience  that  good  SiO^  films  have  negligible  electron  trapping 
states  so  that  electron  currents  leave  behind  negligible  stored 
negative  charge. 

G.  The  P-N  Junction-Diode  Substrate-Detector  Experiments 

One  of  the  goals  of  our  research  has  been  to  identify 
the  sign  of  the  current-carrier  in  the  oxide  during  breakdown. 

This  goal  has  been  achieved  through  the  use  of  a  shallow-junction 
P-N  diode  for  the  substrate,  instead  of  a  homogeneous  wafer,  as 
illustrated  in  Fig.  4.  In  addition  to  being  a  low-leakage 
Junction,  the  main  requirement  on  the  junction  is  that  it  be 
shallow,  i.e.,  sufficiently  close  to  the  interface  that  a 
minority  carrier,  injected  into  the  diode  from  the  oxide,  will 
diffuse  to  the  junction  with  negligible  probability  of 
recombination  before  reaching  the  junction.  In  effect,  the  total 
structure  may  be  viewed  as  a  transistor,  with  the  oxide  film 
serving  as  the  emitter  and  the  diode-proper  furnishing  the  base 
and  collector.  In  the  polarity  used  for.  -F.ig.  Mb)  (negative 
corona) ,  if  the  dominant  current  carrier  in  the  oxide  during 
breakdown  is  holes ,  tunneling  in  from  the  p-region  of  the  diode , 
then  the  only  current  observed  in  the  external  circuitry  will 
be  that  in  the  circuit  connected  to  the  p-region  (the  current 
labelled  I  in  Fig.  Mb}).  Conversely,  in  terms  of  the  transistor 


SHIELD 


SHALLOW  JUNCTION 
SUBSTRATE 


a)  SHALLOW -JUNCTION  DIODE  EXPERIMENT 

DIODE  AREA  :  4X7  mm  ,  REVERSE  DARK  CURRENT  «  |x|0“9amp 
JUNCTION  DEPTH :  5000 A;  DOPING  DENSITY  LO'Mo^cm"3 
^substrate  :  1 "  lOG-cm ;  SIO2  THICKNESS:  IOOOA  OVER 
DIODE ,  6000  A  ELSEWHERE. 


Vs  (volt*) 

b)  NEGATIVE -CORONA  RESULTS  :  Vi, 

SHALLOW-JUNCTION  DIODE  EXPERIMENT ;  RESULTS  FOR  NEGATIVE  CORONA 

FIGURE  4 


-  19 


analog  cited  above,  if  the  dominant  current  in  the  oxide  were 

electrons ,  tunneling  in  from  the  negatively  charged  states  on 

the  exposed  surface  of  the  oxide,  then  one  might  expect  the 

only  external  current  observed  to  be  that  flowing  in  the  circuit 

connected  to  the  n-region  (the  current  labelled  I  in  Eig.  4(b)), 

except  for  a  very  much  smaller  hole-circuit  current  I  due  to 

P 

a  small  amount  of  recombination  in  the  p-region.  The  actual 
measured  currents,  presented  in  the  Fig.,  correspond  to  neither 
situation.  In  e^iect,  sizable  currents  are  observed  in  both 
circuits  with  1^  -  -2  1^.  This  result  is  readily  explained 
in  terms  of  the  known  energy-band  structure  of  the  Si-SiO^ 
interface , as  presented  in  Fig.  1.  Electrons  entering  the 
silicon  from  the  oxide  are  initially  very  hot;  in  fact,  these 
electrons  lie  at  least  3.1  ev  above  the  conduction  band 
minimum  in  the  silicon  -  more,  if  the  electrons  are  also 
hot  in  the  oxide.  Thus,  the  entering  electrons  have  an  excess 
energy  at  least  three  times  the  silicon  bandgap .  The  mean  free 
path  for  creation  of  an  electron-hole  pair  on  the  silicon  by 
such  a  hot  electron  is  exceedingly  short,  probably  less  than 
50A.  (Note  that  the  emission  of  one  or  two  optical,  or  acoustical, 
phonons  by  the  hot  electron,  prior  to  the  impact-ionization  event, 
is  of  no  consequence).  ThuB ,  with  near-certainty  the  one  hot 
electi'on  quickly  becomes  two  electrons  and  one  hole.  The  two 
electrons  are  collected  by  the  n-region,  after  diffusing  to  the 
Junction,  and  the  hole  is  ohmically  relaxed  out  of  the  p-region. 

The  current  data  in  the  Fig.  quite  accurately  reflect  this  picture. 

Starting  from  3  ev,  the  first-generation  progeny  (the  first 
electron-hole  pair)  both  have  excess  kinetic  energy.  However, 
on  the  average  the  parent  electron  and  the  two  progeny  will 
evenly  divide  up  the  total  available  excess  of  2  ev  leaving  each 
of  the  three  without  sufficient  energy  for  another  impact- 
ionization  event.  Note  that  the  rather  accurate  two-to-one 
ratio  of  ll^j/l^  is  rather  cogent  evidence  against  hot  electrons 
in  the  oxide,  despite  the  high  breakdown  field  -  1.3  x  IQ1  V/cm. 

I 
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If  the  electrons  in  the  oxide  had  been  heated  to  2  or  3  ev 
before  entering  the  silicon,  the  corresponding  excess  energy 
in  the  oxide  would  have  been  4  or  5  ev  (subtracting  out  the 
silicon  bandgap)  instead  of  2  ev  and  subsequent  electron-hole 
progeny  would  become  almost  a  certainty,  leading  to  I  /  |l  | 
values  closer  to  unity  than  the  0.5  actually  observed,  and  to 
Ip/ ^ox  =  *p^  (  1^!-  instead  of  the  ~1  actually  observed. 

We  see  that  the  shallow-junction  diode  experiment  not 
only  established  that,  for  negative  corona  in  dry  air,  the 
breakdown  carriers  in  the  oxido  are  tunneling  electrons  from  the 
surface,  but  it  has  also  furnished  very  strong  evidence  that  the 
electrons  in  the  oxide  are  not  strongly  heated  up  by  the  very 

(  7 

high  field,  1  x  10  V/cm.  The  direct  demonstration  of  hole- 
creation  in  the  silicon  near  the  interface  has  important  bearing 
cn  an  interesting  phenomenon  to  be  discussed  below,  and  that  is 
the  build-up  of  positive  charge  in  the  oxide  above  a  certain 
current  level  (~  5  x  10  ^  amp/cm^)  with  negative  corona.  It 

will  be  argued  below  that  this  build-up  is  probably  due  to  the 
tunneling  of  the  impact-generated  hot  holes  from  the  Si  into 
the  oxide. 

With  positive  corona  in  dry  air  the  diode  experiment 
yields  only  one  circuit  current,  that  in  the  circuit  going  to 
the  base  region  of  the  analog  transistor,  i.e.,  the  n-region 
of  the  diode.  From  this  it  is  concluded  that  the  breakdown 
carriers  in  the  oxide  are  electrons  which  have  tunneled  in  from 
the  silicon.  The  C-V  measurements  after  charging  yield  no 
flat-band  shift,  indicating  negligible  trapped  charge  in  the 
oxide— -  completely  consistent  with  the  above  interpretation.  (A 
dominant  hole  current  in  SiO^  will  always  yield  trapped  positive 
charge).  Further  evidence  confirming  electron  current  injection 
comes  from  a  detailed  study  of  the  decaying  surface  potential 
after  charging,  as  discussed  below.  These  measurements,  which 
are  extremely  accurate,  are  in  excellent  agreement  with  a 
theoretical  model  based  on  Fowler-Nordheim  tunneling  of  electrons 
from  the  conduction  band  of  Si  into  the  SiOg. 
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In  addition  to  the  above  purposes  which  the  Junction- 
diodes  have  so  usefully  served,  they  can  also  be  very  useful  in 
obtaining  other  kinds  of  valuable  information.  Virtually  any 
source  of  electron-hole  excitation  in  the  silicon  can  be 
detected  by  the  diodes.  For  example,  we  were  able  to  demonstrate , 
using  the  diodes,  that  light  emitted  in  the  region  of  the  corona 
discharge  around  the  needle  point  and  reaching  the  exposed 
insulator  surface  had  a  negligible  effect  on  our  measurements. 

(in  this  demonstration  the  shield  was  biased  to  cut  off  the  ion 
current  to  the  oxide). 

H.  Preliminary  Results  with  Dry  Nitrogen 

The  results  discussed  up  to  this  point  were  obtained 
in  the  environmental  chamber  in  a  dry  air  ambient.  Some 
preliminary  experiments  were  carried  out  in  a  dry  N2  ambient  in 
a  bell  Jar.  With  positive  corona  the  breakdown  field  in  the 

i 

oxide  was  the  same  as  in  the  dry  air  experiment.  However, 
with  negative  corona  the  breakdown  field  differed  substantially 
from  the  dry-air  value.  These  results  are  pre-eminently  sensible 
in  view  of  our  interpretation  of  the  dry-air  experiments.  So 
long  as  breakdown  proceeds  by  tunnel-injection  of  electrons  from 
the  silicon  into  the  oxide  (positive  corona)  the  breakdown  field 
of  the  oxide  will  obviously  be  independent  of  the  discharge 
ambient.  On  the  other  hand,  if  breakdown  proceeds  from  tunneling 
of  electrons  out  of  surface  states  (adsorbed  ions  from  the 
negative  corona) ,  it  will  be  sensitive  to  the  nature  of  the 
adsorbed  ion,  which  may  be  expected  to  vary  as  the  ambient  is 
changed.  A  systematic  exploration  of  the  breakdown  behavior  of 
the  oxide  in  a  variety  of  gaseous  ambients  is  expected  to  shed 
more  light  on  these  phenomena. 

I .  Steady-State  I-V  Characteristics 

Current-voltage  characteristics  for  the  non-destructive 
breakdown  of  the  oxide  film  induced  by  a  corona  (both  +  and  -  ) 
in  dry  air  are  exhibited  in  Fig.  5.  These  characteristics  were 
obtained  under  steady-state  charging  conditions  using  the 
comparison  technique  previously  described.  The  lower  levels  of 
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BREAKDOWN  OF  Si02  ON  Si  ( POSITIVE  AND  NEGATIVE  CORONA 
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current  are  controlled  by  a  grid,  as  seen  in  Fig.  2,  and  the 
higher  levels  are  controlled  by  varying  the  voltage  on  the 
corona  needle.  The  relatively  large  error  bars  at  the  lower 
levels  of  current  are  a  manifestation  of  the  slow  response  of 
the  insulator,  as  previously  discussed.  At  higher  current 
levels  the  error  bars  are  narrower,  but  still  there  is  not 
sufficient  resolution  to  give  truly  quantitative  results,  i.e. 
results  sufficiently  precise  to  test  a  theoretical  model  for 
I  vs  V.  However,  the  resolution  is  sufficient  to  define  a 
'breakdown  field'.  Further,  there  is  a  substantial  overlap  of 
this  region  with  the  same  region  of  I  vs  V  as  measured  for  the 
decaying  surface  potential  by  the  more  precise  Kelvin  probe 
method,  as  previously  discussed.  The  two  methods  give  the  same 
results . 

Note  that  for  a  given  sign  of  corona  the  same  I  vs  V 
curve  is  obtained  for  both  n-  and  p-silicon  substrates.  For 
the  negative  corona  this  is  an  expected  result  since  the 
breakdown  is  initiated  at  the  exposed,  charged  surface  of  the 
oxide  (by  tunneling  of  electrons  out  of  the  adsorbed  ions  into 
the  oxide  conduction  band).  For  the  positive  corona  the 
coincidence  of  I  vs  V  curves  for  both  n-  and  p-substrates 
conveys  an  interesting  piece  of  information.  The  n-substrate 
is  in  strong  accumulation  and  the  p-substrate  in  deep  inversion. 
Our  results  show  that,  at  room  temperature,  the  pair-generation 
rate  in  the  p-Si  (near  or  at  the  interface)  is  sufficient  to 
supply  the  oxide  current  without  widening  of  the  depletion  layer 
Just  behind  the  inversion  region,  i.e.  without  driving  the 
p-silicon  into  deep  depletion.  By  lowering  the  substrate 
temperature,  the  pair-generation  rate  can  be  made  sufficiently 
slow  (  it  is  a  thermally  activated  process  which  decreases 
exponentially  with  temperature  as  exp  (-EQ/kT))  that  the 
p-silicon  must  go  into  deep  depletion.  At  sufficiently  low 
temperature  only  steady-state  avalanching  in  the  deep-depletion 
region  can  supply  the  oxide  current.  These  effects  have 
already  been  observed,  at  room  temperature,  in  analogous 
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experiments  on  MNOS  structures  (with  very  thin  oxides) 

We  plan  to  look  for  them  in  corona  experiments  on  pure  SiO 
films  at  lower  substrate  temperatures. 

J •  Thickness  Dependence  of  the  Breakdown  Field 

SiOg  film  breakdown  fields  have  been  obtained,  through 

I  vs  V  curves  as  in  Fig.  5,  for  both  +  and  -  corona  in  dry  air, 

over  a  range  of  oxide  thicknesses  between  approximately  500 

and  3000A.  The  data  points,  plotted  in  Fig.  6,  show  a 

negligible  thickness  dependence  over  this  thickness  range.  (The 

o 

apparent  small  drop  at  2700A  is  attributed  to  a  larger 
systematic  error  at  higher  voltages).  This  observed  thickness 
independence  is  important  for  the  theoretical  modeling  of  the 
breakdown-  Rather  similar  results  (thickness  independence)  have 
been  obtained  for  AlgO^  films  on  Si  (the  A1£0  being  deposited 
by  pyrolysis  ),  with  Au  electrodes,  in  an  MOS  capacitor  configura¬ 
tion,  by  Carnes  and  Duffy^,  els  seen  in  Fig.  7  (a),  reproduced  from 
Reference  5,  Shown  in  Fig.  7(b)6,  reproduced  ^rom  Reference  6,  is 
a  collection  of  breakdown  data  for  AlgO^  films  taken  from  a 
number  of  investigators.  The  most  obvious  thing  about  this 
Fig.  is  the  state  of  confusion  as  regards  the  experimental 
breakdown  situation  for  Al^  films.  Most  (possibly  all)  of 
this  confusion  must  derive  from  the  wide  variety  of  deposition 
methods  used  in  the  preparation  of  the  films. 

K.  Electron  Tunneling  Induced  by  Positive  Corona 

Careful  measurements,  using  the  vibrating  Kelvin  probe, 
were  made  of  the  decaying  surface  potential  of  the  oxide , 
following  charging  by  a  positive  corona.  The  decay  was  followed 
over  a  large  time  span,  ranging  from  a  few  seconds  to  many 
hours,  as  shown  in  Fig.  8.  The  positive  corona  polarity  is  of 
particular  interest  because  in  this  polarity  the  geometry  of 
the  experiment  lends  itself  to  detailed  theoretical  modeling, 
namely  on  the  basis  of  Fowler-Nordheim  tunneling  from  the 
conduction  band  of  the  silicon  into  the  oxide.  (In  the  negative 
corona  polarity,  as  we  have  already  seen  from  the  shallow- 
junction  diode  experiments,  we  are  dealing  with  electrons 
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tunneling  out  of  surface  states  (adsorbed  ions)  into  the  oxide. 
Since  a  detailed  knowledge  of  these  surface  states  is  lacking 
at  present  -  we  know  nothing  of  the  trapped-electron  wave 
function,  cr  even  its  binding  energy  -  a  detailed  theory  is  not 
possible  for  the  present.  For  this  reason  the  decay  after 
negative  corona  charging  is  not  discussed  here).  In  order  to 
interpret  the  positive-corona  decay  curve  exhibited  in  Fig.  8 
it  is  necessary  first  to  sketch  the  background  theory. 

First  we  note  that,  since  there  is  no  observed  trapping 
cf  electrons  m  the  oxide,  the  field  across  the  oxide  is  quite 
uniform,  and  accurately  given  by  its  nominal  value:  £  =  V  /d, 
Since  V  is  a  known,  measured  function  of  time,  £  is  also  a 
known  function  of  time:  C,  -  fc(t). 

The  relations  between  the  conduction-current  density 
J  through  the  oxide,  the  surface  charge  Q  per  unit  area  on  the 
oxide,  and  the  oxide  field  £  are 


J(t)  *  -  ^  Q(t)  =  -  e  £  Gt) 


(1) 


where  c  is  the  static  dielectric  constant  of  the  oxide,  and  we 
have  used  Gauss'  theorem. 

For  J  we  use  the  Fowler-Nordheim  tunneling  expression, 
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here  omitting  the  correction  terms: 


J  ■  b£  2  exp(-a/£  ) 


(2) 


with 


0  =  \  * 3/2  “d  6  ■  isfsq-  sr-  * (3J 


with  $  the  barrier  height,  as  defined  in  Fig.  9,  mQX  the 
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electron  effective  mass  in  the  oxide  conduction  hand,  m_.  the 
properly  averaged  (orient  at  ion- dependent)  effective  mass  in 
the  silicon  conduction  band,  "fc  =  Planck's  constant/27r  ,  and  q 
the  charge  on  the  electron. 

Substitution  of  (2)  into  (l)  yields  the  differential 

equation: 


exp(  -  ct /£  )  =  -  e  §- 


(4) 


or 


a 

exp  z 

E? 

which  integrates 


dL=  -  dt 

directly  to 


(5) 


4  =  -  £n  A(t  +  t  )  ,  A  =  —  .  (6) 

£  a  o  e 

The  constant  of  integration  tQ  is  determined  from  (6)  evaluated 
at  t  =  o  : 


~  =  -  In  At 

tQ  “  0 


(7) 


Equations  (6)  and  (7)  yield  the  oxide  field  t  as  a  function  of 
time  t.  Note  that  for  t  >>  t  the  predicted  plot  of  l/£  vs  £n(t} 
is  asymptotically  linear,  as  observed  experimentally  (see  Fig. 9). 
Further,  the  slope  of  the  asymptotic  straight  line  gives 
precisely  the  constant  a  appearing  in  the  Fowler-Nordheim 
expression  (2).  Finally,  the  experimental  plot  yields  a  value 
for  A,  from  which  the  Fowler-Nordheim  pre-exponential  constant  3 
is  also  determined. 

In  processing  the  experimental  data  points  (circles) 
of  Fig.  9»  first  the  linear  portion  is  extrapolated  down  to 
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t  =  0,  at  which  time  the  measured  field  yields  the  constant 

t  ,  using  (7).  When  tQ  is  added  to  each  value  of  t  the  resultant 

plot  (crosses)  in  Fig.  9  is  an  accurate  straight  line  extending 

over  four  decades  in  time  (corresponding  roughly  to  four  decades 

in  current),  and  following  the  predicted  form,  Eq.  (6).  The 

measured  slope  of  this  line  gives  a,  which,  according  to  (3) 

3/2 

is  related  to  m  <p  .  Using  for  <f>  its  recently  measured 
8 

value  we  finally  obtain: 

m 

(  <f>  =  3.1  ev)  —  =  0.39  •  (8) 

m 


Note  that  this  value  for  m  is  in  good  agreement  with 

ox  o 

previously  published  results  of  Lenzlinger  and  Snow  ,  who 

studied  tunneling  in  MOS  capacitor  structures.  The  following 

table  compares  our  results  for  m  with  theirs: 

ox 


m  /m 
ox 

from  <p 

f  Tunneling 

from  Si  : 

0.42 

3.25  ev 

Lenzlinger 

"  A1  : 

0.39 

3.2 

and  Snow 

L  " 

"  Mg  : 

0.48 

2.4 

Our  result 

it 

"  Si  : 

0.35 

3.25* 

(*  Used 

for  purpose 

of  comparison 

only) 

A  theoretical  value  for  A,  the  constant  appearing  in 
Eq.  (6), is  obtained  from  a  more  refined  analysis  of  Fowler- 
Nordheim  tunneling  which  will  be  presented  in  a  more  detailed 
report.  Here  we  state  only  the  result,  namely 


theoretical 


=  1.3  x  10 


15 


sec 


-1 


(9) 


Experimental  values  for  A,  obtained  for  three  different  oxide 
thicknesses  grown  on  n-Si ,  are: 


d  (A)  : 

A  (sec  ^)  ; 


587 

1.8  x  1015 


1000 

2.5  x  10ll+ 


2580 

1.4  x  1015 
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Considering  the  simplicity  of  the  model  this  is  quite 
satisfactory  agreement  with  (9).  It  is  to  be  recalled  that  the 
Fowler-Uordheim  analysis  is  worked  out  for  tunneling  from  a 
highly-degenerate  Sommerfeld  (non-interacting)  electron  gas  into 
the  vacuum.  In  the  present  set  of  experiments  there  is  strong 
band-bending  of  the  silicon  at  the  interface.  Thus  the  degenerate 
electro  gas  at  the  interface  is  actually  held  in  a  tri¬ 
angular  well  narrow  enough  to  produce  significant  space-quantization 
effects.  The  latter  have  been  studied  extensively  by  Stern10 
and  others  at  the  IBM  Research  Laboratories ,  without  considering 
tunneling.  We  have  extended  this  analytical  work  to  include 
tunneling  into  the  oxide,  as  outlined  in  Section  8  of  this 
report . 

In  concluding  this  discussion  of  decay  experiments 
we  would  like  to  emphasize  the  accuracy  and  reliability  of  the 
technique.  The  Kelvin  probe  measurements  have  an  accuracy  of 
+  0.1%.  A  potential  is  directly  measured  rather  than  a  current; 
thus,  weak  spots  in  the  oxide,  of  small  area,  have  a  negligible 
influence  on  the  overall  result.  This  is  in  marked  contrast 
to  MOS  capacitor  measurements  where  current  is  measured  and 
where,  consequently,  weak  spots  can  lead  to  large  errors.  Likewise, 
edge  effects  and  potential  spreading  associated  with  metallic 
contacts  and  high  fields  are  avoided  with  the  present 
technique.  Measurements  on  different  samples  cut  from  the 
same  wafer  give  results  that  agree  to  within  0.1%.  Different 
thicknesses  of  oxide  have  the  same  slope  to  within  1%.  Our 
conclusion  is  that  the  cofona-  charging,  potential-decay  technique 
promises  to  be  a  powerful  technique  for  further  studies  of 
tunneling. 

L'  Positive  Charging  Induced  by  Negative  Corona 

A  final  subject  remains  to  be  discussed  in  this  section, 
and  that  is  the  observation  of  positive  charging  of  the  oxide 
under  negative  corona,  as  noted  on  Fig.  5.  A  typical  charging 
rate,  plotted  as  flat-band  voltage  shift  vs  time  is  shown  in 
Fig.  10.  The  amount  of  positive  charge  corresponding  to 


VFB  =  "20V  ani  d  =  is  Q  '  5  x  1012cnf2  (assuming  the 

charge  is  trapped  near  the  interface,  as  is  certainly  reasonable 
in  the  present  problem.  If  the  charge  were  trapped  near  the 
exposed  sunace,  the  corresponding  Q  would  be  higher/.  The 
simplest  interpretation  of  this  charging  phenomenon  is  that  it 
is  due  to  the  tunneling  of  holes  from  the  Si  into  the  Si02  at 
the  very  high  fields,  exceeding  I  x  107  V/cm.  We  proceed  to 
analyze  the  charging  on  this  basis,  the  model  being  schematically 
illustrated  in  Fig.  11. 

The  fiexd  in  the  oxide  at  the  interface,  Bj  = 
can  be  written,  using  Gauss'  Theorem,  as 


Q't;  being  the  charge  trapped  at  (or  near)  the  interface  at  time 
t,  and  being  the  equivalent  chmic  field,  ^  -  V/d.  (V  is 
constant  in  time,  since  these  are  steady-state  corona  experiments). 

In  the  spirit  of  the  tunneling  model,  we  represent  the 
hole  current  Jfa  entering  the  oxide  from  the  silicon  by  a  Fowler- 
Nordheim-like  expression: 


T  .  dfl(t)  „  ®h 

'  «  *  c  exp  -  ,  (u> 

neglecting  the  pre-exponentiai  field-dependence,  i.e.,  C  #  Ci£./, 
since  the  model  is  not  precise  enough  (at  this  stage)  to 
warrant  such  a  complication.  Note  that  in  writing  ^11)  we  have 
assumed  that  all  injected  holes  are  trapped  in  the  oxide 
(presumably  near  the  interface;.  Again,  departure  from  this 
assumption  would  bring  in  unnecessary  complications  at  this  stage 
of  the  game. 

Substituting  for  £j(t)  from  (10)  into  (ll)  gives  the 
differential  equation  for  Q: 


SCHEMATIC  ENERGY- LAND  DIAGRAM 
FOR  TUNNELING  AND  TRAPPING  OF  HOLES  IN  Si02 


FIGURE  II 


(12) 


a& 

dt 


C  exp  - 


Since ,  even  at  the  upper  end  of  the  V^g  vs  t  plot  in 
Fig.  10,  the  reduction  in  the  field  at  the  interface  is  only 
about  15!  there  is  no  significant  loss  of  accuracy  in  expanding 
the  exponent  in  (12),  giving  the  new  equ  -on 


=  C  exp  -(  'lt1-  (l  +  -S—  )}  =  J  exp  -  j-  £  ,  (13) 

dt  fc0  o 

where  J  =  C  exp  (-<*,/£  )  is  the  initial  hole  flux, 
o  no 

Equation  113)  is  readily  integrated  to  give 

Q(t)  =  in  +  1)  >  (l^) 

with 


where  the  boundary  condition:  Q(o)  =  o,  has  been  used. 

The  corresponding  flat-band  voltage  shift  is 

b  2d 

v  „  _  QUK.  >  -  — —  in  (-  +  1)  (16) 

FB  e  t 

In  view  of  this  theoretical  expression,  the  exper¬ 
imental  data  of  Fig.  10  were  re-plotted  on  a  semi -log  plot 
exhibited  in  Fig.  12.  It  is  seen  that  quite  good  linearity  holds 
over  the  full  two  decades  in  time  for  which  the  data  were 
available.  From  the  slope  of  this  line,  using  (l6) ,  a  value  for 
is  obtained,  namely: 
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p  ^ 

q,  =  (3.7  +  0.8)  x  10  V/cm  •  (IT) 

h 

The  corresponding  Fowler-No rdheim  constant  for  electrons, 
obtained  from  the  surface-potential  decay  measurements,  Fig.  8, 
is : 


a  =  2-3  x  108  V/cm  .  (18) 

e 

From  (17)  and  (l8)  we  see  that,  at  most,  ~  2- 

The  theoretical  value  for  this  same  ratio,  assuminfc 

tunneling  of  thermal  holes  (for  &h)  and  thermal  electrons 

( for  a  )  is : 
e 


theor . 


(^oxjh  jl/2 

m 

ox,e 


hv  3/2 


e 


21 


mox,h 'l/2 
m 


ox,e 


(19) 


where  we  have  used  the  values  -  ^-.8  ev  and  -  3.1  ev  (see 
Fig ,  l).  We  have  already  seen  that  m  -  0.4  m.  Although 

D  OX  jS 

m  .is  not  known,  a  reasonable  guess  is  that  m  , /m  ■'  1,  as 
ox,h  ox»n 

would  follow  if  holes  in  the  oxide  are  either  strongly  self- 

trapped  (a  speculative  concept)  or  drift  in  an  equivalently  very 

narrow  band  (very  high  density  of  states  at  the  top  of  the  valence 

band).  Finally,  note  that  the  calculated  Fowler-No rdheim 

Q  O 

J  -  10  A/cm  is  quite  consistent  with  the  measured  sample 

o  7  2  . 

current,  J  -  5  x  10-1  A/cm  already  determined  to  be  carried 

by  electrons,  whereas  the  calculated  Fowler-Nordheim  t  is 

about  five  orders  of  magnitude  larger  than  the  experimental  x, 

gotten  from  Fig.  12.  All  of  these  observations  suggest  strongly 

that  tunneling  of  hot  holes  generated  in  the  Si  by  the 

incoming  electrons  is  responsible  for  the  charging  of  the  oxide. 

This  hypothesis  will  be  incisively  tested  experimentally. 
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3*  Locally  Destructive  ,  Self-Quenching  Breakdown  in  MIS 
Capacitor  Structures  (D.  Y.  Yang  collaborating) 

A.  Summary 

Work  in  continuing  on  the  self-quenching  breakdown  (SQB) 
studies  of  pyrolytically  deposit6d  Al^  films  on  Si,  and  such 
studies  have  now  also  been  initiated  on  thermally  grown  SiO^  films 
on  Si.  Voltage-time  studies  of  breakdown  events  on  the  Al^O^ 
structure,  with  a  Au  electrode,  consistently  give  an  order-of- 
magnitude  longer  time  scale  for  breakdown  with  the  substrate  in 
inversion  (hundreds  of  nanoseconds),  than  in  accumulation  (tens 
of  nanoseconds).  Charge-storage  (C-V)  measurements  on  the 
Al^O^  structures  suggest  that ,  whatever  the  polarity  of  the  Au 
gate,  the  breakdown  is  initiated  at  the  Si  interface.  Photo¬ 
micrographs  of  the  Au  gate  show  marked  differences  in  the 
topography  of  the  localized  breakdown  areas  between  the 
accumulation  and  inversion  situations.  The  breakdown  voltage  vs 
time  behavior  of  SiO^  structures  follows  a  pattern  very 
similar  to  that  of  the  Al^O^  structures,  with  an  order-of- magnitude 
longer  time  scale  for  inversion  as  compared  to  accumulation. 

C-V  measurements  indicate  negligible  charge  storage  induced  by 
negative  gate  voltage  and  either  large-scale  positive  charge 
storage  and/or  production  of  new  interface  states  induced  by 
positive  gate  voltage.  SEM  studies  of  the  Si02  structures 
after  breakdown  have  shown  marked  differences  ir  the  topography 
of  the  breakdown  spots  produced  by  positive  gate  polarity  as 
compared  to  negative  gate  polarity. 

B.  The  Breakdown  Voltage  vs  Time  Behavior  for  AlgO^ 


Structures 

The  pyrolytic  Alo0  films  have  been  studied  over  a  thick- 

ness  range  300  to  2000A  ,  on  both  p-Si  and  n-Si  substrates,  of 

resistivity  approximately  1  fi-cm,  and  with  gold  metallization 

0  2 

of  thickness  =  1000A  and  area  =0.01  cm  .  A  constant  voltage 
was  applied  to  the  MAS  (metal-A^O^-  Si)  structure  through  a 
series,  current-limiting  resistor  1^.  Following  the  initiation 
of  breakdown,  the  voltage  across  the  film  was  monitored  on  a 
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Tektronix  51*6  oscilloscope.  The  breakdown  pulse  was  used  to 
trigger  the  oscilloscope  sweep.  (However,  because  of  the 
built-in  oscilloscope  delay  line  a  portion  of  the  signal  before 
triggering  can  still  be  seen).  Typical  voltage-time  character¬ 
istics  are  exhibited  in  Figs.  13,  ll*  and  15  which  are  direct 
tracing  of  photographs  of  the  oscilloscope  screen.  Upon 
breakdown,  the  voltage  deceases  sharply,  then  climbs  back 
again,  relatively  slowly,  (on  a  time  scale  determined  by  the 
external  circuitry),  to  the  full  applied  voltage.  The  most 
straight-forward  interpretation  of  the  voltage  drop  is  that 
it  is  produced  by  discharge  of  the  electrostatic  energy  stored 
in  the  capacitor  through  the  breakdown  spot.  The  turn-around 
point  in  the  V  vs  t  trace  presumably  corresponds  to  the  instant 
of  time  at  which  either  evaporation,  or  melting  flow,  or  blow¬ 
out  of  the  gate  in  the  immediate  neightborhood  of  the  breakdown 
volume  produces  an  open  circuit  between  the  intact  area  of  the 
gate  and  the  breakdown  volume.  It  is  this  breakdown-produced 
open  circuit  that  we  call  ' self-quenching' . (It  has  previously 
been  called  'self-healing'.  From  the  gross,  overall  capacitance 
point  of  view,  this  is  appropriate  terminology.  From  the 
microscopic  point  of  view,  there  is  severe  local  damage,  and 
'self-quenching'  seems  a  more  appropriate  description).  The 
recovery  of  the  voltage  across  the  capacitor  is  determined  by  the 
RC  time  constant,  with  P  the  current-limiting  series  resistance 
and  C  the  MAS  capacitance. 

Figures  13(a)  and  (b)  show  the  V  vs  t  characteristic 
o 

for  a  film  thickness  ~  300A,  on  a  p-Si  substrate.  In  (a)  the 
Si  is  in  accumulation  (gate  negative);  in  (b)  the  Si  is  in 
inversion  (gate  positive).  Figure  13(a)  is  a  single-shot  break¬ 
down;  Fig.  13(b)  is  a  multiple  exposure  exhibiting  several  break¬ 
downs.  As  expected,  there  is  some  fluctuation  in  behavior  from 
breakdown  to  breakdown.  Comparing  (a)  and  (b),  it  is  seen  that 
the  time  scale  for  the  breakdown  is  a  full  order  of  magnitude 
longer  with  the  substrate  in  inversion  than  in  accumulation,  an 
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effect  first  reported  by  Klein  in  196611  in  studies  of  SiOg  films 
on  Si.  Concomitantly,  the  fractional  voltage  drop  across  the 
oxide  is  much  larger  for  accumulation  than  it  is  for  inversion. 
Comparable  results  are  obtained  with  a  thicker  Al^O^  film, 
d  =  2000A  ,  on  p-Si,  Fig.  l4,  where  (a)  is  taken  in  accumulation 
and  (b)  and  (c)  in  inversion,  (b)  with  light  illumination  and  (c) 
in  the  dark.  We  note  further,  that  with  the  substrate  in 
inversion  a  microscope  lamp  shining  on  the  sample  through  the 
gate,  Fig.  14  (b),  can  shorten  the  breakdown  time  as  compared  to 
its  value  in  the  dark,  Fig.  lU(c)  -  an  obviously  reasonable  result 
Note  that  (a)  is  a  single-shot  exposure,  (b)  and  (c)  multiple- 
shot  exposures.  Figure  15  exhibits  the  V  vs  t  transient  for  an 
n-Si  substrate  with  dQx  =  1000A.  Figure  15(a)  shows  the  normal 
fast  breakdown  for  the  substrate  in  accumulation.  Fig.  15(b), 
also  for  accumulation,  is  unusual  in  that  it  exhibits  two  distinct 
time  constants  for  the  recovery  following  self-quenching,  a 
normal,  fast  recovery  initially,  barely  visible  on  the  drawn- 
out  time  scale,  1  psec/div,  followed  by  a  very  slow  recovery 
to  full  voltage.  Figure  15(c)  exhibits  several  breakdown 
transients,  with  the  usual  fluctuations.  Some  characteristic 
numbers  for  the  voltage-time  transients,  derived  from  the 
oscilloscope  photographs,  are  presented  in  Table  I.  (see  page  46). 

If  the  time-constant  t  characterizing  the  voltage  vs 
time  transient  is  plotted  vs  dQx  the  results  obtained  are  as 
shown  in  the  log-log  plot  of  Fig.  l6.  It  is  seen  that  with 
the  substrate  in  accumulation  the  data  is  well  represented  by  the 
simple  relation  t*  l/d^.  This  suggests  that  an  equivalent 
breakdown  resistance  Rg  be  defined  thru 
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Table  I.  Some  Characteristic  Numbers  for  the  Voltage-Time  Transients 


*  #*  ##* 


Xtal  type 

Interface 

d 

ox 

V  , 
appl 

VSQ 

Tsq 

T 

p- 

accum 

0 

300A 

-26V 

-i6v 

-9 

UOxlO  sec 

80xlo“9£ 

p- 

accum 

2000 

-150 

0 

20xl0-9 

20xl0"9 

n- 

accum 

1000 

+95 

+10 

20xl0“9 

—6 

20xl0“9 
, ,  ,  -6 

P- 

inv. 

300 

+30 

+27-5 

1x10 

11x10 

P- 

inv. 

2000 

+150 

+100 

100xl0-9 

400xl0-9 

n- 

inv. 

1000 

-95 

-70 

150xl0“9 

550xl0-9 

* 


V  is  the  voltage  at  the  turn-around  point  in  time 
SQ, 

7  is  the  time  to  turn-around 
SQ 


'•V 

***  t  is  the  extrapolated  time  to  reach  0  voltage,  using  the 
slope  of  the  V  vs  t  curve  at  t  =  0. 


with  A  the  capacitor  area,  and  with  the  subscript  a  added  to 

signify  that  the  substrate  is  in  accumulation.  (Correspondingly, 

the  subscript  i  signifies  that  the  substrate  is  in  inversion). 

From  the  plot  of  Fig.  l6,  Rg  5a  =  100ft.  Using  only  the  two  points 

corresponding  to  d^^  =  1000A  and  ^-ox  ~  2000A,  then  also  in 

inversion,  t  <*  1/d  ,  and  R_  .  *  5Kft.  (There  is  a  point  for 

y  1  OX  Ojl 

d  =  300A  which  is  well  off  the  inversion  straight  line.  How- 
ox 

ever,  it  is  also  to  be  noted  that  thin— oxide  substrates,  in 
inversion,  generally  lead  to  shorts  following  breakdown  and  there¬ 
fore,  in  some  sense,  embody  different  physical  phenomena). 

The  V  vs  t  relation  associated  with  the  simple  RC 
representation  (20)  of  the  discharge  is 


V(t>  ■  exp  - 


V  =  d 
B  ox 


(21) 


where  V  is  the  voltage  across  the  oxide  ,  Vg  the  voltage  at 

breakdown,  Eb  the  thickness-independent  breakdown  field,  and 

R  the  equivalent  breakdown  resistance. 
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The  power  dissipated  in  the  resistance  Rg  is 

,  £  „ 

3v  = 


F  =  lv'dtv=0  If!  V 


.  (dox^B)2  2t 

-  - <—  exp  -{ — 


“b 


D  „  ■  d  \ .  (22) 

FLeA  ox 


Equations  (22)  and  (20)  state  that  the  initial  power  dissipation 
is  proportional  to  the  square  of  the  oxide  thickness,  whereas  the 
time  constant  for  the  dissipation  is  inversely  proportional  to 
the  thickness.  Assuming  that  the  equivalent  resistance  can  be 
associated  with  an  approximately  fixed  volume  of  silicon  (or 
silicon  plus  oxide),  as  suggested  by  the  large  difference  in 
time  scales  between  the  accumulation  and  inversion  situations, 
we  would  infer  that  self-quenching  is  favored  by  thicker  oxides. 
Further,  from  (22)  we  see  that  self-quenching  is  also  favored 
by  accumulation  (smaller  Rg)  as  compared  to  inversion  (larger  Rg). 
Indeed,  we  find  that  with  thin  oxides,  dQx  in  the  range  300  to 
800  A»  and  with  the  substrate  in  inversion,  the  breakdown  almost 
always  leads  to  a  short  rather  than  to  self-quenching.  These 
results  bear  directly  on  Carnes'  observation8,8  that  in  his 
ramp-mode  experiments  it  was  very  difficult  to  obtain  consistent, 
reproducible  behavior  when  the  substrate  was  in  inversion, 
although  the  same  substrate,  ii  accumulation,  always  exhibited 
definable,  measureable  behavior .  With  the  thin  oxides,  and  the 
substrate  in  inversion,  it  was  generally  not  possible  to  obtain 
a  voltage  drop,  following  breakdown,  sufficiently  large  to 
trigger  the  oscilloscope  sweep.  However,  observation  of  the 
sample  through  the  optical  microscope  did  reveal  the  presence 
of  bubble-like  areas  on  the  Au  gate,  but  no  visible  evidence 
of  local  destruction. 

Cl  C-V  Measurements  and  Optical  Photomicrography  of 

AlgO^  Structures 

C-V  measurements  were  made  on  the  Al^  structures  to 
determine  the  charge  storage  in  the  oxide,  both  before  the 
application  of  the  high  field  and  after  the  occurrence  of  self- 


quenching  breakdowns,  Figure  17  exhibits  the  results  for  MAS 

o 

samples  with  dQ^  =  1000A  and  with  an  n-type  Si  substrate.  There 
is  an  initial,  positive  flatband  voltage  of  approximately  15V 
prior  to  the  application  of  high  fields.  With  the  gate  positive, 
the  flatband  voltage  shifted  to  +45V.  On  a  different  capacitor 
on  the  same  sample,  with  the  gate  negative,  the  flatband 
voltage  shifted  in  the  opposite  direction,  from  +15V  to  -25V 
after  breakdown.  There  appears  to  be  no  distortion  in  the 
post-breakdown  C-V  curves.  Other  samples  with  a  p-type  substrate 
generally  exhibited  the  same  features.  A  simple  interpretation 
of  these  results  is  that  breakdown  is  initiated  at  the  Si 
interface  rather  than  at  the  Au  interface.  This  important 
conclusion  remains  to  be  verified  by  independent  experiments. 
Similar  charging  effects  could  be  observed  at  high  fields  below 
the  breakdown  threshold.  Pre-breakdown  charging  of  the  oxide 
might  influence  the  breakdown  via  a  distortion  of  the  applied- 
field  distribution  away  from  the  uniform,  ohmic  distribution. 

Observations  of  self-quenching  breakdown  were  made 
through  an  optical  microscope.  Breakdown  events  occurred 
randomly  in  space  and  were  always  accompanied  by  a  flash  of 
light.  Figure  l8  shows  two  typical  photomicrographs,  (a) 
taken  of  a  p-type  sample  in  accumulation,  and  (b)  taken  of 
an  n-type  sample  in  inversion.  Figure  18(a)  is  typical  of 
breakdowns  with  the  substrate  in  accumulation,  whether  n-or 
p-type,  and  also  for  a  p-type  substrate  in  inversion.  Corr¬ 
esponding  to  each  breakdown  event  there  is  an  outer  circular 
area,  of  diameter  =  10  y  and  an  inner  spot  of  diameter  <  1  p  . 

Figure  18(b)  is  typical  of  breakdowns  with  an  n-type  substrate 
in  inversion.  Corresponding  to  each  breakdown  there  is  an  outer 
area  of  quite  irregular  shape,  with  a  maximum  dimension  on  the 
order  of  20  y  roughly,  or  less.  As  in  Fig.  18(a),  there  is 
also  a  much  smaller  inner  spot .  The  photomicrographs  suggest 
that  there  is  an  inner  core  to  the  discharge  corresponding  to 
a  filament  through  which  the  capacitor  energy  discharges. 
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(b)  n-Si  In  Inversion,  d  =  1000A 

ox 


Fig.  18.  Optical  Photomicrographs  of  Localized  Breakdown  Areas 
on  Au-A^O^-Si  Structures.  (Magnification:  300X) 
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D.  The  Breakdown  Voltage  vs  Time  Behavior  for  SiO„ 
Structures  "  - - — — • 

Breakdown  voltage  vs  time  studies,  similar  to  those 
discussed  in  Section  B  for  Al^  structures,  have  also  been 
made  on  Si02  structures.  The  oxides  are  thermally  grown, 
to  an  approximate  thickness  dQx*  2500A,  on  both  n-  and  p-Si 
substrates  of  resistivity  in  the  range  1-2  ft  cm.  The  metal 
electrode  is  Al,  of  thickness  »  1000A,  and  with  the  area  in 
the  range  10  3  to  10  2  cm2.  Typical  volt age- time  character¬ 
istics  are  exhibited  in  Figs.  19  and  20,  which,  again,  are 
direct  tracings  of  photographs  of  the  oscilloscope  screen. 

Figures  19  (a)  and  (b)  exhibit  the  V  vs  t  character¬ 
istics  for  Si02  on  a  p-Si  substrate.  In  (a),  the  Si  is  in 
accumulation  (gate  negative);  in  (b)  the  Si  is  in  inversion 
(gate  positive).  Figures  20  (a),  (b)  and  (c)  correspond  to 
an  n-type  substrate  -  (a)  to  accumulation.,  (b)  and  (c)  to 
inversion.  The  gate  area  in  (b)  is  double  that  in  (c),  and 
accordingly  the  breakdown  time  constant  is  in  (b)  approximately 
twice  that  in  (-).  As  in  Section  B,  one  can  use  the  measured 
time  constants  to  define  equivalent  breakdown  resistances, 
obtaining  Rg^a  £  75  0  in  accumulation,  and  =  1.5Kft 

in  inversion.  In  this,  as  well  as  all  other  respects,  as 
regards  the  voltage-time  measurements,  the  behavior  of  SiO 
structures  is  very  similar  to  that  of  AlgO  structures.  Thus, 
the  time  scale  for  breakdown  is  about  an  order  of  magnitude 
larger  for  the  inversion  case  than  for  the  accumulation  case, 
and  the  fractional  voltage  drop  across  the  oxide  is  much 
larger  in  the  accumulation  case,  90-  100*,  than  for  the  inversion 
case,  k0-60%  (with  large  fluctuations). 

E.  C-V  Measurements  on  SiC>2  Structures 

Charge  storage  in  the  oxide  was  monitored  by  C-V 
measurements  as  exhibited  in  Fig.  21  for  a  p-type  substrate, 
and  Fig.  22  for  an  n-type  substrate.  The  solid  lines  denote 
the  initial  C-V  areas,  measured  before  the  samples  were  exposed 
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to  high  fields.  THey  indicate  that  there  is  negligible  charge 
storage  prior  to  the  application  of  the  high  fields.  In  Fig.  21(a), 
the  inversion  case, it  is  seen  (the  dashed  curve)  that  positive 
gate  bias,  inducing  breakdown,  produce  a  large  negative  flatband 
shift  of  about  TO  volts,  corresponding  to  a  volume  density  of 
positive  charge  (near  the  Si  interface)  10  /cm  .  A  substantial- 
distortion  of  the  C-V  curve  is  produced,  accompanied  by  large 
hysteresis.  These  striking  results  are  consistent  with  field- 
induced  tunneling  of  holes  from  the  A1  and  their  subsequent 
transport  through  the  oxide  and  trapping  near  the  interface. 

(However,  based  on  the  corona  studies  reported  in  Section  2L, 
we  see  that  this  observation  does  not,  by  itself,  establish  the 
sign  of  the  dominant  carrier  in  the  oxide  during  breakdown. 

Further  experiments  are  needed  to  settle  this  point).  Note 
that  the  magnitude  of  the  capacitance  is  somewhat  reduced  due 
to  loss  of  gate  area  through  repeated  breakdown  events-  The 
magnitude  of  the  nominal  breakdown  field  is  8  x  10  V/cm.  In 
Fig.  21(b),  the  accumulation  case,  it  is  seen  (the  dashed  curve) 
that  negative  gate  bias  leading  to  breakdown  results  in  no  charge 
storage  whatever.  This  strongly  suggests  breakdown  initiated  at 
the  A1  gate  with  the  injected  electrons  as  the  dominant  carriers. 

For  both  n-  and  p- substrates ,  whether  in  accumulation  or  inversion, 
there  is  a  steady  occurrence  of  breakdown  events  at  the^field  strength 
8  x  10^  V/cm,  and  an  onset  of  such  events  at  6,6-7.0x10  V/cm. 

Figures  22(a)  and  (b)  exhibit  comparable  plots  for 
an  n-type  substrate,  (a)  corresponding  to  the  gate  positive 
(substrate  in  accumulation)  and  (b)  to  the  gate  negative  (substrate 
in  inversion).  Although  the  gate-positive  case  does  not  produce 
a  large  flatband  shift  as  in  Fig.  21(a),  there  is  very  substantial 
distortion  with  some  hysteresis, evident  in  the  post -breakdown 
((dashed)  curves  of  Fig.  22(a),  indicating  a  probably  large 
production  of  interface  states.  Figure  22(b)  strongly  resembles 
Fig.  21  (b)  in  that  the  post -breakdown  (dashed)  curve  is  very 
little  changed  from  the  pre-breakdown  (solid)  curve,  indicating 
negligible  charge  storage  and  negligible  surface-state  production 
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due  to  the  breakdown. 

Studies  are  also  under  way  of  pre-breakdown  current 

flow  and  charge  storage  at  high  fields.  With  the  gate  of  a 

p-Si  sample  held  positive  at  a  nominal  field  of  7-4  x  10  V/cm 

6  2 

the  current  decayed  from  an  initial  value  of  10  A/ cm  by  one 
order  of  magnitude  over  a  24  hour  period.  Concomitantly,  the 
C-V  plots  exhibited  a  slow  drift  towards  negative  voltage  - 
a  10  volt  shift  after  1.5  hours  and  a  30  volt  shift  after 
24  hours.  These  pre-breakdown  effects  are  currently  being 
studied  systematically. 

F.  SEM  Topographic  Studies  of  SiO^  Structures 

Topographical  studies  of  the  local  breakdown  areas  on 

SiO^  structures  were  made  in  the  Princeton  SEM,  with  results  for 

p-type  samples  exhibited  in  Figs.  23  and  24.  A  notable  feature 

which  the  SiO^  breakdown  spots  have  in  common  with  the  Al^O^ 

breakdown  spots,  Fig.l8,  is  the  existence  of  an  inner  core 

surrounded  by  a  much  larger  area  from  whicn  the  gate  material  has 

been  removed.  We  have  already  noted  the  strong  polarity  dependence 

of  the  behavior  of  the  C-V  plots.  Figs.  21(a)  and  22(a)-  as  compared 

to  Figs.  21(b)  and  22(b).  Polarity  dependence  also  shows  up  in 

the  SEM  photographs.  With  the  gate  positive,  Figs.  23(a)  -  (c), 

the  edge  of  the  breakdown  area  is  quite  round  and  very  smooth 

in  appearance,  as  if  molten  during  the  course  of  the  breakdown. 

Figure  23(a)  shows  a  few  overlapping  spots  and  23(b)  a  single 

spot.  Figure  23(c)  exhibits  the  structure  after  stch-off  of  the 

gate.  The  inner  core  of  the  spot  appears  vivir  y  exposed  as  a 

hole  'burned  through'  the  oxide.  For  all  three  photographs  of 

-2  2 

Fig.  23  the  total  gate  area  was  the  same,  namely  1.3  x  10  cm  . 
Figures  24  (a)  and  (b)  exhibit  breakdown  spots  for  the  opposite 
polarity,  gate  negative,  and  for  two  different  gate  areas, 

1.3  x  10-2cm2  in  (a)  and  1.3  x  10_3cm2  in  (b).  In  both  cases, 
the  breakdown  areas  have  a  quite  different  appearance  form  those 
seen  in  Fig.  23  (gate  positive).  Those  in  Fig.  24  are  smaller 
and  have  much  rougher  edges,  as  if  gate  material  were  either  blown 


(a)  Scan  of  Gate 
(5800X) 


(b)  Scan  of  Gate 
(5800X) 


SEM  Photographs  of  Localized  Breakdown  Areas  on  Al-SiO 
Structures;  p-Si  Substrate,  Gate  Positive. (Gate  Area  =' 
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(b)  Gate  Area  =  1.3  x  10  cm(58,000X) 


Fig.  24.  SEM  Photographs  of  Localized  Breakdown  Areas  on  Gate  of 
Al-Si02~Si  Structures;  p-Si  Substrate,  Gate  Negative. 
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out  or  exceedingly  rapidly  quenched  from  a  molten  state.  The 
spot  in  (a)  is  much  larger  than  either  of  the  two  in  (b), 
corresponding  to  the  order-of-magnitude  larger  gate  area  in 
(a),  and  correspondingly  larger  stored  capacitor  energy.  The 
SEM  photographs  of  Figs.  23  and  24  were  taken  on  p-type  substrates. 
With  n-type  substrates,  negative  gate  voltage  (substrate  in 
inversion)  produces  breakdown  spots  of  exceedingly  irregular 
shape  as  already  noted  in  the  optical  photograph  of  Fig.  18(b). 

With  the  oxide  completely  etched  off,  penetration  of  the  inner 

core  of  the  breakdown  into  a  p-Si  substrate  has  been  seen  in 
the  SEM. 

An  interesting  phenomenon  was  discovered  accidently 
when  four  dots  on  p-type  substrates  were  inadvertently 
scratched.  With  positive  gate  polarity,  the  breakdowns  were 
clearly  initiated  along  the  line  of  a  scratch,  whereas  with  fi¬ 

ne  gat  ive  gate  polarity  the  scratch  appear  to  have  no  effect 
on  the  spatially  random  pattern  of  breakdown  spots.  These 
same  results  were  also  found  to  hold  for  n-type  samples. 

After  etch-off  of  the  oxide,  no  sign  of  the  scratches  could 
be  seen  on  the  substrate. 

k'  Techniques  for  Probing  the  Solid-Solid  Tn+.prf,na 

A*  Background  Discussion 

When  an  electron  beam  strikes  a  solid  surface  there  are 
a  number  of  "signal  channels"  activated  which,  if  properly 
collected  and  analyzed,  can  offer  useful  analytical  information 
about  the  sample.  These  are  illustrated  schematically  in  Fig.  25. 
The  electron  beam  loses  energy  when  passing  through  a  sample  by 
a  number  of  mechanisms.  These  include  plasmon  generation,  phonon 
generation,  "knock  on"  damage,  and  ionization  processes, both 
localized  and  non-localized  in  nature.  These  energy  losses 
limit  the  range  of  the  electron  beam  and,  at  the  same  time,  give 
rise  to  analytical  information  resulting  from  the  relaxation  of 
the  absorbed  energy.  In  most  of  the  existing  electron-optical 
installations,  provisions  are  made  to  gather  the  information  from 
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only  one  or  two  of  these  signal  channels.  In  order  to  deal  in 
the  most  efficient  way  with  a  particular  signal,  specialized 
instruments  have  been  developed  such  as  the  Auger  electron 
spectrometer,  the  electron  microprobe,  the  transmission  electron 
micrc scope  and  the  scanning  elect  - on  microscope.  The  most 
complete  specification  of  a  material  can  be  obtained  only  by 
collecting  the  signal  from  a  maximum  number  of  activated 
signal  channels  in  one  instrument.  An  electron  beam  technique 
has  been  developed  which  can  probe  structure  at  the  interface 
between  a  thin  insulating  film  and  a  semiconductor  substrate 
utilizing  the  simultaneous  collection  of  signals  from  more 
than  one  signal  channel. 

B.  Differential  Detection-Volume  Spectroscopy 

The  information  contained  in  the  various  signal 
channels  provides  analytical  data  for  the  volume  of  material 
in  which  a  given  signal  is  generated.  This  volume  we  will 
refer  to  as  the  "detection  volume".  It  i~  schematically 
illustrated  in  Fig.  26  for  a  number  of  signal  channels  being 
excited  by  the  same  incident  electron  beam.  The  exact  sizes 
of  the  detection  volumes  for  different  signal  channels 
depend  both  on  the  nature  of  the  specimen  under  examination 
and  on  the  beam  voltage,  which  together  determine  the  range 
of  the  primary  electrons,  and  the  angle  of  incidence  of  the 
incoming  electron  beam,  which  determines  the  surface  normal 
component  of  that  range.  The  depth  to  which  this  detection 
volume  extends  can  therefore  be  altered,  by  changing  the 
pertinent  parameters,  from  a  very  shallow  penetration  on  the 
order  of  a  few  angstroms  to  a  penetration  of  several  microns. 
The  range  R  of  the  electrons  in  solids  can  be  expressed  by 
the  following  approximate  expression: 


where  k  and  n  are  constants  for  a  given  material  and  Eq  is  the 
primary  beam  energy.  This  expression  is  accurate  to  within 
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•*  10  per  cent  in  the  range  of  incident  energies  between  2  and 
50  kev  for  metallic  absorbers .  It  is  less  accurate  for 
materials  of  lower  electron  density.  Since  the  values  of  n 
and  k  are  not  readily  available,  and  the  range  of  electrons 
m  thin  films  of  varying  density  is  not  expected  to  follow 
this  equation  with  great  accuracy,  an  alternate  method  of 
estimating  the  range  can  be  used.  Assuming  that  the  density 
of  a  material  is  proportional  to  the  volume  electron  density 
we  may  determine  a  range-density  product  for  a  given  incident 
beam  energy  and  then  estimate  the  range  for  any  other  material 
by  dividing  the  product  by  the  appropriate  density.  Estimates 
of  primary  electron  range  using  this  range-density  product  can 
be  made  with  an  accuracy  of  +_  30% .  The  errors  introduced  by 
the  change  from  one  material  to  another  are  smaller  than  those 
due  to  the  changes  in  the  densities  of  materials  when  we  are 
examining  thin  films.  In  order  to  examine  the  interface  between 
a  semiconductor  substrate  and  a  thin  insulating  film  we  can 
alter  the  range  of  the  electrons  (and  thereby  the  detection 
volume)  so  that  the  volume  sweeps  through  the  interface  region. 
When  the  detection  volume  includes  part  of  the  interface  and 
substrate  regions  we  may  record  differences  in  the  signal 
channels  and  attribute  them  to  the  differences  in  the  material 
encompassed  by  the  detection  volume,  i.e.,  the  interface  and 
the  substrate.  The  information  gathered  in  this  way  is  of 
course  an  integral  of  the  particular  signal  channel  over  the 
entire  detection  volume  so  that  the  sensitivity,  resolution, 
etc.  of  the  technique  are  best  when  the  insulating  film  on  the 
substrate  is  as  thin  as  possible  and  there  is  little,  if  any, 
metallization  on  the  sample  surface.  For  an  analysis  of  the 
interface  between  an  insulating  film  and  its  substrate  the 
following  measurements  are  made:  1.  The  chemical  composition 
of  the  oxide  is  determined  using  X-ray  fluorescence  with  the 
sample  at  a  high  tilt  angle.  The  incident  beam  voltage  must  be 
sufficiently  high  that  the  energy  levels  of  interest  will  be 
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excited  by  the  beam, and  the  penetration  depth  is  controlled 
by  using  a  high  angle  cl  incidence  sc  that  the  electron  range 
normal  to  the  surface  is  less  than  the  thickness  of  the  oxide. 

2.  The  chemical  composition  of  the  interface  region  and  the 
substrate  material  is  determined  using  X-ray  fluorescence  ,  as 
above,  by  decreasing  the  angle  of  tilt  so  that  the  detection 
volume  includes  the  interface  and  tart  cf  the  substrate  region. 
'The  spectra  are  then  normalized  to  those  obtained  from  the 
insulating  film  alone  Urigh  tilt  angle;  to  determine  that 
portion  of  the  spectra  which  is  due  to  the  material  at  the 
interface  and  in  the  substrate.  A  further  increase  in  the 
penetration  will  allow  some  separation  of  the  spectra  due 
to  the  insulating  film  and  the  interface  region  from  that  due 
to  the  substrate.  An  example  of  such  a  measurement  is  shown  m 
Fig.  27.  While  these  measurements  are  not  yet  quantitative 
in  depth  resolution,  they  ofie/  information  about  the  critical 
region  near  the  interface  which  is  net  a  ,  suable  from  any 
other  nondestructive  technique  at  'his  time  Our  expectation 
is  that  future  deveiopmei  t  of  this  technique  will  yield 
semiquantitative  or  quantitative  results  3.  The  structure 
of  the  oxide  can  be  probed  by  using  an  induced-current  mode 
of  imaging  where  the  insulating  film  is  sandwiched  between 
the  substrate  and  a  thir.  metai  layer.  The  secondary  emission 
of  the  sample  at  low  voltages  may  charge  the  specimen  either 
positively  or  negatively  to  provide  a  sell -biased  induced- 
current  image.  As  the  incident  angle,  and  electron  energy, 
are  changed  we  can  form  the  images  by  the  transport  of  charge 
due  to  the  eiectron-beam-mduced  carriers  at  varying  depth  in 
the  specimen.  An  example  of  an  image  formed  m  this  way  is 
shown  in  Fig.  28.  Figure  28(a;  is  an  image  formed  with  incident 
electron  beam  energy  of  30  Kev  and  an  incident  angle  of  30 
degrees  using  the  secondary  electrons  to  modulate  the  intensity. 
Figure  28(b)  is  a  similar  image  formed  with  the  same  beam  energy 
and  incident  angle  but  the  induced  current  signal  has  been  used 
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Figure  28(a)  Secondary  electron  image  showing  surface  structure. 


Figure  28(b)  Film  and  interface  region  imaged  with  induced  current 
signal.  Structure  not  visible  in  the  surface  image 
(a)  above  is  noted. 


Figure  28.  SEM  secondary  electron  and  induced  current  images  of 

A1203  on  Si.  The  beam  energy  is  30  KV  and  the  incident 
angle  is  30°. 
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to  modulate  the  image  intensity.  The  electron  beam  penetration 
into  the  oxide,  after  compensation  for  the  stopping  power  of  the 
top  electrode,  is  on  the  order  of  1300  angstroms  and  the  impact 
area  of  the  beam  includes  the  interface  region.  The  structure 
visible  in  the  induced  current  image  can  be  attributed  to  the 
oxide  or  interface  region  since  it  is  not  visible  in  the  secondary 
electron  image  which  has  information  about  the  surface  layer 
only.  Figure  29  shows  a  lower  magnification  of  the  same  area. 

In  Figure  29(a),  a  secondary  electron  image,  taken  with  a  beam 
energy  of  15  Kev,  the  striations  on  the  surface  of  the  sample 
are  visible  as  in  Figure  28.  Figure  29(8)  shows  that  at  lower 
accelerating  voltages  the  lateral  nonuniformities  in  the  charge 
stored  at  the  interface  and  in  the  oxide  region  can  be  imaged 
using  the  induced  current  signal.  Spatial  resolution  on  the  ord  r 

n 

of  200  A  is  obtainable  with  this  technique.  Although  interpre¬ 
tation  of  these  images  is  not  yet  well  understood,  they  offer 
a  method  of  observing  the  structure  at  the  interface  which  has 
not  yet  been  exploited,  and  work  is  currently  underway  to 
quantify  this  method  for  the  study  of  thin  insulating  layers 
on  semiconducting  substrates. 

5 .  Charge-Discharge  Studies  of  Electron-Trappling  in  AH- Implanted 
SiO^  (Noble  Johnson  collaborating) 

A .  Background  Discussion  and  Summary 

It  is  a  well-known  property  of  clean,  thermally  grown 
Si02  films  on  Si  that  they  are  relatively  free  of  electron  traps, 

In  consequence,  electrons  are  easily  removed  from  the  oxide  at 
high  electric  fields,  there  generally  being  no  impediment  to 
their  egress  at  the  oxide  interface  with  either  the  silicon 
substrate  or  the  gate.  On  the  other  hand,  holes  are  very 
strongly  trapped  in  the  oxide,  either  throughout  the  bulk  or  at 

# 

Work  supported  in  part,  also,  by  ONR  under  Contract  NOOO14-67-A-O151-0O35 


V 


Figure  29(b)  Film  and  interface  region  imaged  with  induced  current 

signal.  Two-dimensional  charge  distribution  is  visible. 


Figure  29.  3 EM  secondary  electron  and  induced  current  images  of 

Al^O^  on  Si.  The  beam  energy  is  14  Kev  and  the  incident 

angle  is  40°. 
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the  SiOg-Si  interface.  The  very  weak  trapping  of  electrons  and 
the  very  strong  trapping  of  holes  is  almost  certainly  the  origin  of 
the  radiation  sensitivity  of  SiOg  films  on  Si.  One  logical 
approach  to  radiation-hardening  of  the  SiOg  is  controlled  doping 
of  the  oxide  in  such  a  manner  as  to  produce  large  numbers  of 
electron  traps  and  hence  strong  trapping  of  electrons  in  the  oxide. 

A  convenient  means  of  controlled  doping  of  the  oxide  is  ion 
implantation.  Since  aluminum  implantation  has  shown  promise  for 
radiation  hardening  of  the  oxide,  we  have  been  studying  the 
electron-trapping  behavior  of  Al-implanted  SiOg  films  on  silicon 
by  charge-discharge  techniques.  These  studies  are  not  concerned 
directly  with  radiation-hardness,  but  rather  with  the  fundamental 
electronic  behavior  of  SiOg  as  modified  by  the  implantation- 
dcping  with  Al.  A  brief,  preliminary  report  of  these  studies  was 
given  in  the  First  Semiannual  Technical  Report.  A  more  extensive 
repor*  is  presented  here. 

The  experiments  have  conclusively  demonstrated  that  unannealed 
aluminum-implanted  SiOg  possesses  electron  traps  not  found  in  the 
undoped  material.  Optical  excitation  of  trapped  electrons  indicates 
that  the  traps  are  located  in  the  lower  half  of  the  forbidden 
energy  gap  of  SiOg.  A  model  of  electron  trapping  has  been  form¬ 
ulated  which  illustrates  the  dynamic  features  of  the  trapping 
process  and  correctly  predicts  the  time  dependence  of  the  flat- 
band  voltage  and  photocurrent  during  electron  photoinjection. 
Application  of  this  model  to  the  experimental  results  indicates 
that  essentially  all  of  the  photoinjected  electrons  are  trapped 
in  the  oxide.  An  analytical  solution  of  the  rate  equations  has 
been  obtained  by  suitable  approximations. 

B .  Experimental  Procedure  and  Results 

The  samples  used  were  MOS  structures  in  which  the  oxide 

o 

layer  was  grown  to  a  thickness  of  1400  A  in  dry  oxygen  on  the 

(lOO)  face  of  a  3  ft-cm  n-type  silicon  wafer.  The  oxide  of  the 

test  samples  was  implanted  with  20  keV  aluminum  ions  to  a  fluence 
14  -2 

of  10  cm  .  A  control  sample  was  not  implanted.  The  profile 
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of  the  implanted  aluminum  ions  has  been  published  elsewhere. 


Semi-transparent  gold  field  plates  were  vacuum  deposited  on  the 
exposed  oxide  surface.  The  implanted  samples,  as  received 
(unannealed) ,  showed  a  flat-hand  voltage  typically  of  -15  volts 
(see  curve  labelled  1  in  Fig.  5  of  the  First  Semiannual  Report). 

An  initial  exposure  of  the  implanted  samples  to  UV  light 
with  the  contacts  short-circuited  reduced  the  flat-band  voltage 
approximately  to  zero.  After  this  optical  annealing,  a  positive 
voltage  was  applied  to  the  gate  and  the  sample  was  exposed  to 
b.77  eV  photons  from  a  high- intensity  monochromator.  This 
exposure  induced  an  internal  photo-emi3sion  of  electrons  from  the 
silicon  into  the  oxide.  Up  to  the  highest  applied  voltage,  35 
volts,  no  steady-state  current  could  be  observed,  only  a  transient- 
current  which  always  decayed  to  zero.  Following  the  current 
decay,  C-V  measurements,  made  at  1  MHz,  showed  a  positive  flat- 
band  voltage,  indicating  negative  charge  storage,  with  the  magnitude 
of  the  flat-band  voltage  being  essentially  equal  to  the  applied 
gate  voltage.  The  flat-band  voltage  could  then  be  reduced 
approximately  to  zero  by  either  optical  or  thermal  annealing  with 
the  contacts  short-circuited.  Optical  annealing  required  photon 
energies  greater  than  about  b.6  eV.  Thermal  annealing  could  be 
accomplished  at  temperatures  below  350°C.  This  experiment  was 
repeated,  with  reproducible  results,  over  the  entire  range  of 
gate  voltages . 

In  marked  contrast  with  the  foregoing  results  obtained  with 
Al-implanted  samples,  the  non-implanted  control  sample,  when 
subjected  to  the  same  experiment,  showed  no  detectable  storage 
of  negative  charge  and  the  photo-induced  current  showed  no  signs 
of  decaying  over  a  long  period  of  time. 

The  dynamics  of  the  negative  charging  of  the  Al-implanted 

samples  were  investigated  by  shutting  the  light  off,  and  so 

interrupting  the  photoinjection  at  intervals  during  the  charging 

process,  and  measuring  the  C-V  curve.  The  resulting  time 

evolution  of  this  curve  for  a  gate  voltage  of  10  volts  is  shown 

in  Fig.  30.  In  Fig.  31(a)  is  shown  a  plot  of  V  -  V_  vs.  time, 

g  FB 

where  V  is  the  gate  voltage  and  V  is  the  flat-band  voltage, 

6  i*  ij 


C  (pF) 


e  During  Electron  Photoinjection  with  Vq  =  IO  volts 
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The  corresponding  photocurrent,  as  measured  in  the  external 
circuit,  is  shown  in  Fig.  3l(b) .  Similar  results  were  obtained 
for  all  other  gate  voltages  used  in  the  study. 

A  straightforward  interpretation  of  the  above— described 
results  is  as  follows:  The  aluminum  implantation  into  the  SiO^ 
created  deep  electron  traps  in  the  insulator.  Electrons 
photoinjected  into  the  insulator  were  captured  into  these 
traps,  building  up  a  negative  space  charge  which  gradually  reduced 
the  electric  field  at  the  Si-SiOg  interface  to  zero,  concomitantly 
reducing  the  injected  photocurrent  to  zero.  The  trapped  negative 
charge  shifts  the  flat— band  voltage  in  the  positive  direction, 
and  a  shift  in  flat-band  voltage  eq_ual  to  the  applied  gate 
voltage  corresponds  to  the  reduction  in  electric  field  at  the 
interface  to  zero.  In  the  next  section  we  present  a  model  of 
electron  trapping  which  correctly  predicts  the  time  dependence 
of  VFB  and  of  the  Phot°c^rent  measured  in  the  external  circuit 
during  electron  photoinjection.  The  shapes  of  the  curves  in 
Figs.  31(a)  and  (b)  reflect  the  dependence  of  the  photoinjected 

current  on  the  intensity  of  the  electric  field  at  the  Si-SiO 
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interface.  A  tentative  conclusion  resulting  from  application 

of  the  model  to  the  data  of  Figs.  31(a)  and  (b)  is  that  essentially 

all  of  the  photoinjected  electrons  are  trapped  in  the  oxide,  the 

measured  gate  current  being  almost  entirely  caused  by  displacement 

current  in  the  oxide.  In  this  connection  it  should  be  noted 

that  at  the  largest  gate  voltage  used  in  our  experiments,  namely 

35  volts,  the  number  of  trapped  charges  was  approximately 
13  -2 

10  cm  in  the  steady  state,  whereas  the  number  of  implanted 
aluminum  atoms  was  an  order  of  magnitude  greater  than  this . 

The  indicated  position  of  the  centroid  of  the  space-charge 

O 

distribution  of  trapped  charge  was  approximately  620  A  from  the 
metal-oxide  interface. 

C .  Analysis  and  Numerical  Results 

In  this  section  we  examine  two  aspects  of  the  exper¬ 
imental  results:  first,  we  analyze  the  relationships  between  the 


i 
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steady-state  magnitudes  of  photoemitted,  trapped,  and  externally 
measured  charge,  and  shew  by  comparison  with  the  known  aluminum 
impurity  profile  that  it  is  likely  that  all  of  the  photoemitted 
electrons  were  trapped  in  the  oxide;  then  we  analyze  the  time 
evolution  of  the  flat-band  voltage  and  the  externally  measured 
current  during  the  trap-charging  transient, 
i )  Steady-State  Analysis 

The  total  charge  photoinjected  into  the  oxide  at  the 
Si-SiO^  interface,  ,  can  he  expressed  as  the  sum  of  two 

components:  the  charge  trapped  in  the  oxide,  Q  ,  and  the  charge 

O  X 

that  is  drafted  completely  through  the  oxide  and  into  the 

metallic  gate,  Q,,  .  Thus: 

thru 


Q.  =  Q 

u.nj  t)x 


+  *th 


ru 


(21+) 


The  charge  that  has  been  transported  through  the  external  circuit, 
on  the  other  hand,  consists  of  the  sum  of  the  conduction-current 
component  Q  ru  and  a  displacement-current  component  which  is 
equal  to  weighted  by  the  factor  1  -  x/L,  where  x  is  the  dis¬ 
tance  from  the  Metal-SiO^  interface  to  the  centroid  of  the  steady- 
state  distribution  of  trapped  charge  and  L  is  the  thickness  of 
the  oxide ,  i .e . : 


Q  =  (l  -  £)Q  +  CL, 

ext  L  ox  thru 


(25) 


Finally,  the  change  in  flatband  voltage  caused  by  Q  is  given  by 


Z\V. 


FB 


(26) 


where  f  is  the  absolute  dielectric  constant  of  the  oxide. 

By  elimination  of  the  quantity  Q  between  Eqs .  (25)  and 
(26)  we  obtain  the  following  expression  for  the  ratio  of  the 
through-going  change  to  the  charge  measured  in  the  external 
circuit: 
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thru 

ext 


1  - 


eAV. 


FB 


QextL 
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(27) 


By  elimination  of  Qqx  and  Q^.  ^  between  Eqs.  (24),  (25), 
and  (26),  we  obtain  an  expression  for  the  ratio  of  the  charge 
trapped  in  the  oxide  to  the  total  photoinjected  charge: 


Q 

Q 


ox 

inj 


_ L/x _ 

«extL/E  4VFB  +  1 


(28) 


We  now  apply  the  foregoing  relationships  to  the  data  of 

the  present  experiment.  By  integration  of  the  data  for  current 

—7  p 

given  •' n  Fig.  31(b)  we  find  Q  =  2.9  x  10  coulomb/cm  .  Using 
L  =  1.4  x  10  ^  cm,  AV  =  10  volts,  and  a  relative  dielectric 

PD 

constant  of  3.8  for  the  oxide,  we  can  plot  the  ratios  given  oy 
Eqs,  (27)  and  (28)  as  functions  of  the  ratio  x/L.  The  results 
are  shown  in  Fig.  32,  together  with  a  graph  of  the  aluminum 
impurity  profile.12 

Inasmuch  as  Q  /Q.  cannot  exceed  unity,  and  cannot 

ox  inj  J  thru 

be  negative,  the  results  plotted  in  Fig.  32  show  that  the  centroid 

of  steady-state  charge  distribution  is  restricted  to  the  range 

0,45  ^  x/L  j_  1.  Furthermore,  the  tail  of  the  aluminum  impurity 

profile  to  the  right  of  x/L  =0.45  contains  only  about  one-tenth 

of  the  total  number  of  implanted  A1  atoms.  If  we  assume  that  the 

location  of  the  centroid  is  somewhnt  to  the  right  of  x/L  =  0.45, 

say  at  x/L  =  0.6,  we  must  suppose  that  the  implantation-related 

13  -2 

traps  caused  by  fewer  than  10  cm  impurity  atoms  were  able  to 

capture  three-fourths  of  the  electrons  (Q  /Q.  ,  ^  0.75)  while 

ox  inj 

the  traps  on  the  left,  which  were  the  result  of  an  order-of- 
magnitude  greater  number  of  impurity  atoms,  were  unable  to 
capture  the  remaining  25  percent  of  the  carriers.  This  argument 
makes  it  appear  highly  likely  that  essentially  every  emitted 
electron  -was  captured,  and  that  the  location  of  the  centroid  of 
charge  distribution  in  the  final  steady  state  was  at  x/L  =  0.45. 


5£'5"  L _  .  Si  02  (L  =  1400 A)  Silicon 

I  Al -  Implanted ,  20  KeV,IOl4cm2  *"j  Substrate 

Fig.  32.  Plots  of  Relative  Trapped  Charge  and  Relative 
Transmitted  Charge  vs  the  Relative  Trapped-Charge 
^entroid  Position.  Also  Plotted  is  the  Implanted 
Aluminum  Profile  (Ref.  12  ). 
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li )  Transient  Analysis 

We  next  analyze  the  dynamics  of  the  trap-charging 

process,  assuming  that  the  location  of  the  centroid  of  trapped 

charge  distribution  is  not  a  function  of  the  time.  This  assumption 

is  made  purely  for  mathematical  convenience;  its  relationship 

to  the  ongoing  physics  is  discussed  below.  When  a  positive 

voltage  is  impressed  on  the  gate  of  the  MOS  structure  and 

electrons  are  photcemitted  from  the  Si  into  the  SiO^  by  means  of 

optical  excitation,  the  injected  current  is  a  function  of  the 

13 

electric  field  intensity  at  the  Si-SiOg  interface.  The  dependence 
of  phctocurrent  on  electric  field  intensity  for  our  samples, 
using  photons  cf  4,77  eV  energy,  is  shown  in  Fig.  33.  (This 
relationship  was  obtained  by  measurements  made  on  our  unimplanted 
control  sample,  in  which  sample  charging  was  shown  to  be  negligible.) 
As  the  traps  in  the  implanted  sample  capture  electrons,  the  electric 
field  intensity  at  the  interface,  ,  is  gradually  reduced  toward 
zero,  thus  diminishing  the  photoinjec^  -d  current;  meanwhile  the 
trapped  negative  charge  increases  the  flatband  voltage  in  tne 
positive  direction.  The  process  ceases  when  the  electric  field 
in  the  oxide  at  the  interface  is  reduced  to  zero,  at  which  point, 
assuming  negligible  fixed  charges  in  interface  states,  the  bands 
in  the  semiconductor  are  flat;  i.e.,  the  flatband  voltage  has 
been  brought  to  the  value  of  the  gate  voltage. 

Neglecting  the  eflects  of  fixed  charge  at  the  Si-Si02 
interface,  the  flatband  voltage  of  the  MS  structure  can  be 
written  as 

VFB(t)  =  -  7  J  X  °ox(x>t)  331  ’  (29) 

o 

where  0  is  the  metal-semiconductor  work-function  difference, 
ms 

p  is  the  charge  density  in  the  oxide,  and  x  is  measured  from 
ox 

the  metal -insulator  interface.  The  total  charge  per  square 
centimeter,  Qcx>  is 


Dependence  of  the  Electron  -  Current  Photoinjected  from  the  Si 
(£;  is  the  Electric  Field  at  the  Interface) 


-  80  - 


Q  (t  ;  -  '  0  ( x , t )  dx  . 

CX  'OX 

'  o 

The  centroid  of  the  charge  distribution  is  defined  by: 


•L 

!  x  p  (x,t )  dx 
ox 


Ip  (  x,t )  dx 
i  ox 

‘  o 


(30) 


(31) 


or  by  use  c  f  ( 3C ) : 


ox 


P  (x,t) 
ox  ’ 


dx 


(32) 


Using  (321,  we  can  write  (29)  as 


V  ft)  =  0 
FB  ms 


-  —  Q  (t) 

£  OX 


(33) 


Thus,  th’  effect  of  the  distributed  charge  in  the  oxide  is  the 

same  as  would  be  produced  by  a  sheet  of  charge  of  magnitude 

0  legated  at  the  position  of  the  centroid  of  the  distribution, 
ox 

As  was  mentioned  previously,  we  shall  assume  for  simplicity  of 
analysis  that  the  location  of  the  centroid  is  independent  of  the 
time.  Assuming  also  that  all  of  the  photoin jected  charge  is 
trapped  in  th®  oxide,  we  have  the  following  relation: 


dQ 


ox 


dt- 


=  J 


inj 


(3M 


where  J  is  the  density  of  the  photoinjected  current.  This 
inj 

current  iE  a  function  of  the  electric  field  intensity  in  the  oxide 
a+  the  Si-SiOg  interface,  &  ,  in  the  manner  shown  in  Fig.  33. 

For  the  field  intensity  we  can  write 
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where  V  is  the  gate  voltage 
g 


t  (t>  =  i[v  -  0  -  0  +  §•  Q  (t)  1  > 

Csx[r  L  g  ms  s  e  cx 
voltage  and  0  is  the  surface  potential  of 


(  35 ) 


the  semiconductor.  The  current  measured  in  the  external  circuit, 


J  ,  is  given  by 
ext 


Jextft)  =  U  "  ^/L)Jinj(t) 


( ?6) 


An  analytical  solution  can.  be  obtained  by  using  a  piecewise- 
1  inear  approximation  to  the  photoinjected  current  as  a  function  of 
electric  field  intensity,  m  the  manner  shown  in  Fig-  3?.  The 
use  of  two  straight  lines  divides  the  field  dependence  into  a 
high-field  regime  and  a  low-field  regime,  with  the  intersection 
of  the  two  linear  segments  defining  a  transit  ’  field  ^  Hence 
the  photocurrent  density  may  be  written  as 

(  jo  *  ^ 


(High-field  Regime) , 


in  j 


J 


3  A 


k  <C  (Low-field  Regime)  , 

»  Ci  -Cx 

(37) 


,here  J  is  the  intercept  of  the  high-field  linear  segment,  and 
,  and  °  are  the  slopes  of  the  high-  and  low-field  segments, 
respectively  •  For  simplify  the  metal-silicon  work  fnnctron 
difference  and  the  surface  potential  of  the  semiconductor  are 
neglected  in  comparison  with  the  gate  voltage.  From  Eq.  (35) 

P.  depends  only  on  the  applied  gate  bias  and  the  trapped  charge 

2.  _  rm  _ _ X»  _  irn+' 


From  Eq .  ( 33  > ,  depends  only  on  Qc 


Then  for  an  initial 

t  -  — o.  — r -  *  'OX' 

■let -band  voltage  of  zero  the  initial  condition  for  the  total 
rharge  is  Q  (0)  =  0.  Thus,  Eq.  (3k)  can  be  solved  to  yield 
j  (t)  which  when  substituted  into  Eqs .  (35)  and  (3D  yields 
°X  (t) ,  Then  V  (t)  and  J  (t>  are  obtained  from  Eqs  (33) 
andJ(.6),  respectively.  The  following  results  are  then  obtained: 


L 
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In  order  to  make  a  self-consistency  check  with  the  results 
quoted  earlier  in  this  section,  where  the  location  of  the  final 
steady-state  centroid  of  charge  distribution  was  deduced  to  be 
at  x/L  =  0.^5,  we  can  use  Eqs .  (39)  and  (^*0)  to  obtain  the  following 
expression  for  x  (as  obtained  from  our  transient  data): 


.x 
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1  -  (1 


<rV 

g _ 

'  L  J  . 
Vo  ext 
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(UU) 


with 

,  ,  J  L 

J_  *  i_  ( !  t  JE_)  d>5) 

T  **’  Q  V 

Vo  h  h  g 

Substitution  of  the  above-listed  numerical  values  (and  also 

=  li,3  x  102  sec)  into  Eq.  (1*U)  yields  the  result  x/L  =  0.1*9, 

V° 

which  is  only  slightly  in  excess  of  the  value  of  0.U5  which  was 
deduced  from  the  steady  state  analysis  (Fig.  32).  It  is  to  be 
expected  that  the  value  of  x  deduced  from  the  transient  analysis 
should,  if  anything,  be  larger  than  that  obtained  in  the  steady 
state,  for,  as  the  traps  charge,  the  centroid  of  charge  distri¬ 
bution  should  move  away  from  the  Si-SiOg  interface  and  toward 
the  metal  gate;  thus  the  mean  value  of  x  during  he  transient 
should  be  larger  than  the  final  value.  The  discrepancy  between 
the  two  values  of  x  is  small  and  the  difference  is  in  the 
expected  direction;  thus  we  regard  the  resu  s  to  be  satisfactorily 
self-consistent . 

#  ,  .  . 
6 .  Electron-Beam  Injection  and  Damage  Studies  (C.T.  Shih  collaborating  ) 

A.  Background  Discussion 

In  the  initial  studies  with  the  electron  beam, 
summarized  in  the  First  Semiannual  Technical  Report,  we  were 
concerned  primarily  with  the  properties  of  the  beam-induced 
I-V  characteristic  of  the  oxide.,  Typical  results  of  this  phase 
of  the  study  are  presented  in  Figs.  3Ma)  and  (b).  There  is 
If 
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generally  a  temperature-insensitive  square-law-regime,  up  to 

the  point  of  unity  gain  (I  ,  =  I,  ) ,  in  which  furthermore, 

sample  beam 

at  fixed  voltage,  I  I,  .We  interpreted  this  set  of 

sample  beam 

results  in  terms  of  shallow-trapping  of  the  injected  electrons 
in  the  oxide  bulk ,  the  de-trapping  being  due  to  photons 
generated  by  the  beam.  In  order  to  obtain  supporting  evidence 
for  this  interpretation,  extensive  charge-storage  (C-V) 
measurements  were  made,  both  before  and  after  electron-beam 
irradiation.  These  measurements  have  confirmed  negative- 
charge  storage  in  the  oxide  as  a  result  of  the  beam  injection. 
Further,  the  trapping  states  are  created  in  the  oxide  bulk 
by  the  beam,  presumably  within  the  range  of  the  beam. 

Interface  states,  at  the  Si-SiC^  interface,  are  also  created 
by  the  beam  -  undoubtedly  by  photons  generated  by  the  beam  - 
but  these  states  cannot  explain  either  the  I-V  or  C-V  data. 


B .  Charge-Storage  Studies:  C-V  and  Internal  Photoemission 
Measurements 

A  typical  set  of  C-V  curves  is  shown  in  Fig.  35. 

The  capacitance  measurements  are  made  at  high  frequency, 
namely  1  MHz.  The  curve  labelled^ is  taken  on  tne  sample 
prior  to  any  electron  bombardment.  It  indicates  a  small 
initial  positive  charge  in  the  sample.  The  curve  labelled 
@is  measured  following  a  30-minute  exposure,  at  low  temperature 
(107°K)  to  UV  illumination  (hv  =  5  eV)  through  the  gate  with 
the  gate  positive,  still  prior  to  any  electron  bombardment. 

(This  curve  is  actually  independent  of  polarity).  The  very 
small  shift  to  the  right  corresponds  to  a  slight  diminution 
of  the  net  positive  charge,  due  to  trapping  of  the  internally 
photo-injected  electrons.  Electron-trap  densities  in  the 


as-grown  thermal  oxide  generally  do  not  exceed  about 
lU  -3lM5 

10  cm  giving  a  maximum  trapped-electron  sheet- 
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density  below  10  cm  .  This  is  consistent  with  the  tiny 
shift  from^to(^).  The  curve  labelled^)  exhibits  the  C-V 


curve  measured  after  electron -beam  irradiation  for  15  minutes 
at  2.5  kV  with  the  field  in  the  oxide  at  3  x  10°  V/cm.,  with 


Sample  HH8I-I 


ion 


the  gate  negative.  This  C-V  curve  is  obviously  strongly  distorted 
and  also  indicates  large  negative  charge  storage.  Its  general 
shape  is  discussed  further  in  connection  with  Fig.  36.  The 
curve  labelled©  is  a  post-irradiation  curve  taken  after  a 
vacuum  anneal  at  300°C,  for  1  hr,  with  the  gate  negative  at 
hO  V.  A  substantial  recovery  from© toward© is  evident,  both 
in  the  shape  of  the  curve  (much  less  distortion  in0than  in 
0)  and  in  the  amoont  of  shift.  If  the  anneal  had  proceeded 
for  many  hours0would  have  been  driven  back  to©  (this  has  been 
observed) , indicating  the  removal  of  all  negative  charge  trapped 
during  the  electron  irradiation.  After  curve© was  taken,  and 
without  any  further  electron  irradiation  or  annealing,  the 
sample  was  subjected  to  UV  irradiation  (hv  =  5  eV)  for  20  minutes 
at  room  temperature  with  the  results  exhibited  as  curve©  The 
significance  of ©  is  that  it  demonstrates  conclusively  that 
although  the  charge  created  by  the  electron  irradiation  is 
indeed  annealable  at  300°C  (by  thermal  stimulation  out  of  the 
traps)  the  traps  themselves  do  not  appear  to  be  annealable 
( comparing©  with  ©  ),  at  least  at  the  relatively  low  temper¬ 
ature  of  300°C.  The  C-V  measurements  alone  do  not  establish 
the  spatial  location  of  the  electron  traps  in  the  oxide. 

It  is  by  now  well-established  that  almost  any  form  of 
radiation  will  either  directly  or  indirectly  (via  UV  or  higher- 
energy  photons)  generate  interface  states  at  the  Si-SiC>2 
interface.  The  non-penetrating  electron  beam  is  no  exception. 
Using  the  Gray-Brown  technique1^,  involving  C-V  measurements 
over  a  wide  temperature  range  (from  100°  to  300°K,  in  our 
case)  we  have  indeed  shown  that  the  electron  irradiation 
generated  an  order-of-magni'tude  increase  in  the  surface-state 
density  from  2  x  1012  cm-2  to  1  x  1013  cm"2.  However  we  do 
not  see  any  way  in  which  interface  trapping  of  electrons ,  by 
itself,  can  explain  our  I-V  data  summarized  in  Fig.  3*+. 

Further  evidence  of  the  bulk  nature  of  the  electron  trapping 
is  given  below,  in  the  discussion  of  Figs.  ^0  and  4l. 


The  shape  of  the  post-irradiation  C-V  curves Q) and@ in 

Fig.  can  be  understood  by  reference  to  Fig.  36.  In  the 

upper  let''’  inset  we  schematically  illustrate  the  metallization 

dot  with  rwo  areas  indicated,  an  area  exposed  to  the  electron 

beam  'which  is  aper*ured*  and  an  unexposed  area  A.  Purely  for 

simplicity  of  argument  we  assume  A  *  A^-  Further .assume  that 

The  beam  is  uniform  across  A^  (rectangular  profile  in  the 

inset'  Obviously  we  are  dealing  with  two  capacitors  in 

parallel,  areas  A  and  A^^  respectively.  The  unexposed  area  A 

obviously  exhibits  the  same  C-V  curve  before  and  after 

irr .diaT-icn ,  wuich  we  draw  schematically  as  the  dashed  curve 

labelled  C  =  C  in  Fig.  36.  The  exposed  area  A.  exhibits 
b,pre  t> 

a  quite  different  C-V  curve  after  irradiation,  gt  in  Fig. 

■a 6 ,  than  before  irradiation,  C.  .  There  is  a  substantial 
’  b  ,pre 

flat -band  shift  positive  voltages  due  to  stored  negative 

charge  and  seme  distortion  due  to  either  interface-state 

generation  or  lateral  non-uniformity.  The  two  capacitors,  of 

area  A  and  A^  respectively,  are  obviously  in  parallel,  so  that 

the  total  ef f“cti* e  capacitance  is  simply  the  sum  of  the 

separate  -apacitances .  Prior  to  irradiation  C  =C  =  2  C, 

*he  upper  solid  curve  in  Fig.  36,  and  following  irradiation 

C  =  C  *  c  +  r  ,  the  dot -dashed  curve  in  Fig.  36. 

t  t.pcst  b.pcst 

fU;+?  t haT  non-uniformities  in  the  irradiating  beam  will 
inevitably  accentuate  the  distortion  in  the  transition  from 
low  to  high  capacitance  Practically,  the  best  that  can  be 
achieved  with  a  stationary  beam  is  a  Gaussian  profile,  as 
illustrated  in  the  inset  of  Fig.  36,  A  small  lateral  high- 
fiequency  wobble  cf  a  Gaussian  beam  might  then  yield  an 
effectively  uniform  beam.]  Note  the  strong  resemblance  of 
curves  ©and0  in  Fig  35  to  the  schematic  post  curve  in 
Fig.  36- 

The  experimental  C-V  curves  of  Fig.  35  and  schematic 
'C-V  curves  of  Fig  36  are  all  high-frequency  curves.  A 
clearer  separation  of  the  separate  capacitance  effects  of  areas 
A  and  A^  is,  in  principle,  achievable  if  low-frequency  C-V 


Configuration 
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curves  are  measured.  Practically,  it  is  simpler  to  keep  the 

measuring  frequency  fixed  and  rea'-h  the  low-frequency  C-V 

regime  by  either  raising  the  temperature  or  illuminating  the 

Si  cubstrate,  through  the  cxide,with  strongly  absorbed 

^  x8 

( sup’-a-bandgap)  light  *  .  A  comparison  of  the  high-  and 

lew-  frequency  situations  is  given  in  the  schematic  C-V 

diagram  of  Fig.  37 .  The  room-temperature  curves  (T  =  T^) 

labelled  C,  C  and  C  .  are  the  same  as  in  Fig.  36, 

b,pcst  t,post  e  ’ 

era  are  prototype  high-frequency  curves.  The  corresponding 


ru-  es  labelled  C.  ,  C,  .  .  and  ^  ,  are  the 

lew’  b, post, lew  t, post, low 

same  C-V  curves  measured  at  a  sufficiently  higher  temper¬ 
ature  (T  =  T^)  to  represent  prototype  low-frequency  curves. 

The  transition  portion  of  the  high-frequency  C-V  curve  becomes 


a  dip  in  ♦he  ccr responding  low-frequency  C-V  curve.  Further, 


where  a  significant  flat -band  shift  from  C  to  C^  has 

b,post 

occurred  (high-frequency  curves),  as  in  Fig.  36,  the  dips  d^ 
and  d^  in  the  corresponding  1 ow- frequency  curves  are  well- 
sepera* ed ,  as  shown  in  Fig.  37.  The  effects  illustrated 


schematically  in  Fig.  37,  are  seen  experimentally,  as  shown  in 
Fig  38.  Three  post-irradiation  C-V  curves  are  shown, 


corresponding  t0  successively  higher  temperatures  of  measure- 
menT  at  room  temperature  ,  (2)  at  200°C  and  ©at  250°C. 

Bci-h  ©and  ©have  the  high-frequency  appearance  of  curve© 
c:  Fig  35,  or  C^  of  Fig.  36,  whereas  ©exhibits  the 

characteristic  low-frequency  dips  of  Fig.  37.  Curve©  is 
a  post-irradiation,  room-temperature  C-V  curve  following  a 
partial  anneal.  Cuive©is  shifted  leftwards  fron©,  as 
expected.  The  dip  d^,  associated  with  injected  negative 
charge,  is  shifted  leftwards  as  a  result  of  annealing,  the 


shift  being  greater  the  longer  the  annealing  time. 

A  curve  such  as©  in  Fig.  35,  corresponds  to  a  single 
value  of  gate  bias  during  irradiation (in  this  case  — lUO  V). 
Based  upon  the  model  of  Fig.  36,  a  flat -band  voltage  shift 
can  be  associated  with  curve©,  namely  from  that  portion  of 
the  curve  associated  with  '•'he  irradiated  area.  The  flat-band 


.  37.  Effect  of  Temperature  on  C-V  Curves  for  Parallel 
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voltage  shift  so  determined  as  a  function  of  gate  bias  during 
electron  irradiation  (for  a  fixed  dosage  of  3.5  x  10  amp  of  2 
kV  electrons  for  5  hrs  )  is  shown  in  Fig.  39,  icr  irradiation 
at  two  temperatures,  200°K  and  300°K.  The  dip  m  this  curve 
with  increasing  (negative)  gate  voltage  is  not  presently 
understood.  The  large  hysteresis  exhibited  by  the  300°K 
curve  (dashed  portion)  is  quite  normal—  as  I V  !  is  reduced 

o 

in  successive  irradiations,  trapped  electrons  are  not  thermally 

released  and  the  sample  is  not  in  thermal  equilibrium  at  the 

lower  !v  1  A  similar  hysteresis  is  seen  with  the  200°K 
6 

curve  (not  shown) 

The  UV-induced  (hv  =  5-0  eV)  internal  photo-emission 
current  from  the  Si  substrate  into  the  oxide  as  a  function 
of  ga:e  bias  (positive)  is  shown  in  Fig.  ^0,  under  three 
conditions:  curve  ©  is  taken  prior  to  irradiation  by  the 
electron  beam,  curve^^after  irradiation  by  the  beam  and 
curve  ^^after  irradiation,  then  annealing  at  300°C  for  one 
hour.  These  curves  provide  very  strong,  independent  evidence 
for  electron  trapping.  Curve^^has  the  expected  shape 
previously  predicted  and  found  experimentally  by  Powell  and 
Berglund^.  The  post-irradiation  curve  ©is  approximately 
a  translation  ofQby  about  25  V.  This  latter  is  readily 
interpreted  as  the  voltage  required  to  overcome  the  potential 
barrier  created  by  the  trapped  negative  charge  left  m  the 
oxide  bulk  by  the  electron-irradiation.  The  anneal  at  300°C 
for  one  hour  removes  over  half  of  this  trapped  charge  so 
that  curve (3) is  much  closer  to  the  pre-irradiotion  curve  Q 
The  curves  in  Fig.  hO  were  obtained  with  a  slow,  linear 
ramp  voltage  acros-.  the  oxide,  as  indicated  by  the  diagram 
in  the  inset  in  Fig.  *+1.  The  directly  measured  photocurrent 
is  exhibited  as  the  dashed  curves  in  Fig.  Ul.  There  is  a 
hysteresis  associated  with  the  direction  of  the  ramp,  as 
shown.  The  dashed  curve  includes  the  contribution  to  the 
total  current  of  the  constant  displacement  current  associated 
with  the  linear  ramp.  With  the  direction  of  the  ramp  inverted. 


Flat- Bond  Voltage  Shift  vs  Gate  Bias  During  Irradiation 
Different  Temperature 


Photoemission  Current 


Internal  Photoemission 


the  constant  displacement  current  changes  sign,  hence  the 
hysteresis.  The  three  curves  exhibited  in  Fig.  ^0  have  the 
displacement  current  substracted  out,  for  the  sake  of  clarity 
(that  is  why  the  curves  in  Fig.  ^0  all  go  through  the  origin, 
whereas  the  two  dashed  curies  in  Fig  ^1  are  symmetrically 
displaced  from  the  origin). 

More  important  than  the  hysteresis,  which  is  simply  an 
artifact  of  the  measuring  technique,  is  approximate  symmetry 
as  between  positive  and  negative  gate  bias  of  either  dashed 
curve,  after  the  constant  displacement  current  is  substracted 
out.  In  either  polarity  there  is  the  same  magnitude  of 
threshold  voltage  befcre  photo-emitted  current  can  flow  This 
is  precisely  what  one  would  expect  from  a  potential  barrier 
inside  the  oxide  bulk  due  +c  trapped  negative  charge 

Also  exhibited  in  F"g  U1  are  the  quasi-static  C-V 
curves  at  room  temperature  (the  solid  curves),  which,  for  the 


slow  ramp,  are  simply  I  vs 

o  x 


20  21  22 

V  ’  ’  ,  in  the  absence  of  the 

g 

UV.  In  the  upward  direction  of  the  ramp  (-  to  +)  IQy  is  the 

pure  displacement  '’urrent  of  the  ramp  (the  horizontal  line); 

in  the  downward  direction  of  the  ramp  there  is  a  very  rich 

structure  in  I  vs  V  which  is  not  understood  at  present, 
cx  g 

At  elevated  temperatures  such  structure  could  be  attributed 
to  ionic  motions,  but  it  is  difficult  to  invoke  such  ionic 
motions  at  room  temperature 


C.  Suuma ry 

The  outstanding  result  tc-da+e  of  our  electron- 

beam  irradiation  studies  is  the  demonstration  of  stored 

negative  charge  in  the  bulk  following  the  non-penetrating 

irradiation.  Relatively  few  other  studies  have  been  made 

23  2b  25 

utilizing  non -penetrating  electron  beams  ’  ’  and  these 

studies  have  only  reported  positive  charge  storage,  near  the 
Si-SiOg  interface,  following  the  irradiation-  A  major  goal 
of  our  current  program  is  the  reconciliation  of  these 
earlier  results  with  cur  own  results  Our  current  efforts 
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are  concentrated  on  improving  the  uniformity  of  our  electron 
beam  and  designing  means  to  profile  the  team.  We  expect  to 
obtain  *he  greater  uniformity  by  going  over  from  the  standard 
oxide  which  we  have  been  using,  to  a  thoriated  tungsten  cathode. 
Special  phosphor  screens  have  been  fabricated  for  beam¬ 
profiling  purposes,  by  the  Phosphor  Group  at  RCA  Laboratories. 

Theory  oi  Electron-Tunneling  from  Si  into  SiO^  -  Triangular 

Well  Approximation  (Z.  Weinberg  collaborating) 

A .  Background  Discussion 

At  the  high  fields  (>  5  x  10^  V/cm;  see  Fig.  5) 
which  are  required  to  produce  significant  tunneling  currents 
from  Si  into  SiO^  there  is  very  strong  band-bending  in  the 
Si  near  the  interface.  Our  interest  here  is  in  the  polarity 
with  the  Si  negative,  in  which  case  there  is  a  degenerate 
electron  gas  in  the  Si  conduction  band  adjacent  to  the 
interface.  If  the  Si  is  n-type ,  it  is  in  accumulation  at  the 
interface;  if  the  Si  is  p-type ,  it  is  in  inversion.  Whether 
the  Si  is  in  a'-cumulat ion  or  inversion,  the  electron  gas  is 
confined  in  a  very  narrow  well  in  the  direction  perpendicular 
to  the  interface.  This  confinement  leads  to  spatial  quanti¬ 
zation  effects  in  the  perpendicular  direction,  which  have  been 
studied  in  the  inversion  case  by  Stern^,  and  in  the  accumula- 
tion  case  by  Duke.  The  spatial  quantization  effects  have 
been  studied  theoretically  with  various  degrees  of  refinement. 
A  completely  accurate  theoretical  study  requires  the  self- 
consistent  solution  of  the  SchrCdinger  equation  for  all 
electrons  trapped  in  the  well  together  with  the  Poisson 
equation  (since  the  detailed  shape  of  the  well  is  determined 

by  space-charge  considerations).  Such  studies,  which  require 

10 

the  use  of  a  high-speed  digital  computer,  have  been  made. 
However,  in  no  case  which  we  are  aware  of  has  the  spatial - 
quantization  problem  been  looked  at  with  the  aim  of 
calculating  the  tunneling  current  from  the  well  into  the  oxide 
film.  This  is  precisely  the  problem  to  which  we  address 
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ourselves  here.  In  order  to  get  a  reasonable,  analytical 
handle  on  the  problem  we  neglect  the  specific  space  charge 
effects  and  assume  that  the  well  has  a  triangular  shape, 
as  shown  in  Fig.  42.  For  the  case  that  most  of  the  electrons 
m  the  well  occupy  the  lowest  quantized  energy  state,  the 
triangular  well  approximation  is  a  good  one.  Since  this  is 
also  a  case  of  great  practical  interest,  the  results  pre¬ 
sented  here  are  of  direct  use  in  the  analysis  of  tunneling 
currents.  (From  the  diagram  in  Fig.  42  it  is  clear  that  we 
are  also  neglecting  the  well-known  image-force  lowering  of 
the  potential  barrier  by  the  tunneling  electron  This  effect 
produces  only  a  small  change  in  the  tunneling  probability, 
and  can  be  suitably  incorporated  into  the  results  at  a  later 
point . 

B .  The  Mathematical  Calculatips- 

The  ore-electron  (fixed  potential)  quantum- 

mechanical  problem  to  be  solved  is  defined  by  Fig.  42.  The 

electron  in  the  Si  is  trapped  m  the  well  defined  by  the 

potential  step,  height  0  -  3.1  ev  (for  SiOg),  at  the  interface 

and  by  the  bottom  energy  of  the  Si  conduction  band  E  = 

c  ,Si 

-qtx  where  q  is  the  magnitude  of  charge  on  the  electron, 
the  field  in  the  Si  (assumed  constant  near  the  interface)  and 
x  the  position  variable  normal  to  the  interface  which  is 
located  at  x  =  0.  The  barrier  through  which  it  must  tunnel 
is  defined  by  the  potential  step  at  the  interface  and  by  the 
lowest  energy  of  the  SiC>2  conduction  band  Ec  gi0  =  q«5  -  q^x. 

Mathematically  it  is  convenient  to  obtain  the 
solution  in  two  non-overlapping  regions:  Region  I  extending 
from  the  interface  into  the  Si,  and  Region  II  extending  from 

the  interface  into  the  SiO^.  In  each  region  the  solution 

•  .  27 

consists  of  a  suitable  combination  of  Airy  functions  .  The 

separate  solutions  are  then  matched  at  the  interface  to 

provide  current  continuity.  As  one  might  expect  in  a  tunneling 

problem,  the  solution  inside  the  well,  i.e.  Region  I,  is 


I 

I 

I 
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only  slightly  affected  by  the  presence  of  tne  field  in  the 
oxide.  Thus,  the  wave  function  normalisation  can  be  carried 
out  using  only  the  solution  in  Region  I.  This  procedure 
leads  to  a  closed-form  expression  for  the  tunneling  current 
per  electron  out  of  the  ground  state  of  the  well.  The  total 
tunneling  current  mure  include,  additionally,  the  contributions 
from  electrons  in  participating  excited  states  in  the  well. 

Since  the  triangular- well  approximation  is  not  a  good  one 
for  the  study  of  the  excited  levels,  we  can  obtain  only  a 
somewhat  crude  estimate  for  the  tunneling  contributions  from 
the  former  using  Airy-function  solutions. 

We  proceed  to  the  mathematical  discussion: 

Region  I  (Si):  -  K  <  x  <  0,  Field  electron  effective- 
mass  m  . 

The  appropriate  Schrodinger  equation  in  the  normal 
direction  is 

,2.  .  2m  - 

d-  (E  +  xMx)  *  0  .  (46) 

dxd  *  1 

Region  II  (SiC>2):  0  <  x  <  » ,  Field  electron  effective-mass  m2< 
The  appropriate  Schrodinger  equation  in  the  normal 
direction  is 

+  _iL  (E  -  q®  +  q£ x)i|»(x)  =  0.  (47) 

dx  n 

E  in  Eqs.  (46)  and  (47)  is  the  eigen-state  energy  of  the  electron. 

It  is  convenient  to  introduce  dimension-less  variables: 

Region  I:  y1  ■  ad/3  (j^r  +  x)  *  ^  =  •j-'T'  ^ 

Region  II:  yg  =  ad/3  ( +  x) ,  ag  ■  qSg  .  (49) 

Substitution  of  (48)  into  (46),  and  (49)  into 
(47)  reduces  the  two  Eqs.  (46)  and  (47)  to  a  common  form: 
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*  yl,2  *'yI,2)  1  0  '  (50) 

which  is  a  firm  of  Airy’s  Eq  This  equation  has  two  independent 

27 

?oluTicns  which  are  conventionally  written  as  : 

Ai(-y)  and  Bi(-y) 

The  boundary  conditions  appropriate  to  the  present 
tur,r»ling  problem  are  physically  obvious. 

i)  The  wave  function  must  decay  exponentially  going 
toward  -or  in  Region  I,  l  e.  well  beyond  x  =  x^. 

ii )  Ite  va.e  function  must  describe  a  pure  outgoing 
wave  in  Region  II  bey or d  x  *  x^ . 

These  boundary  conditions  enable  us  to  write  the 


appropriate 

solutions  of  (50)  as: 

Regicn  I : 

if  =  a  Ai  (-y1 )  ; 

(51) 

Regi on  I I : 

if  =  T(Bi(-y2)  +  i  Ai(-y2))  . 

(52) 

From  48)  and  (49),  at  x  =  0 

yl. 

1/3  E  .  1/3 

=  a,  -r-  and  y  =  a. 

o  1  2,0  2 

E  -  qgS  _  1/3  i 

“SEf - "2  ^ 

153) 

since  o  *  3-1  ev  E.  With  y  known,  E  is  determined  by 

l  ,o 

(53'.  Using  '  49 )  and  the  known  tunneling  field  on  the  oxide 
(•  6  x  10^  V/cm)  we  obtain  a  good  estimate  for  a^ . 

a^/?  l  2.5  x  109/meter  =  0.25/A  .  (54) 

Using  this  value  in  (53),  taking  j6/£g  =  3.1/6  x  10^ 

1  50A,  we  obtain 
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y?  n  k  10  .  (55) 

It  is  a  property  of  the  Airy  functions  that,  in 
going  from  x2  towards  x  ■  0,  Ai(y2)  decays  exponentially 
whilst  Bi(y2)  increases  exponentially.  Since  y^  >-»  1, 

Ai(  -y^)  can  be  neglected  compared  to  Bi(-y2)  at  y2  =  y^  q  and 
the  matching  of  the  wave  function  at  x  *  0  gives,  using  (51) 
and  (52): 

aAi(-y  )  =  T  Bi(-y0  )  .  (56) 

1,0  2,o 

Matching  of  the  derivative  of  the  wave  function  at 
x  *  0  gives : 


a2/3TB1'(->'2,o) 


(57) 


Because  !y 


2  .o’ 


>>  1  it  is  permissible  to  use  the 


asymptotic  form  of  the  function  Bi(-y2)  in  evaluating  it  at 


y2  =y2,o- 
the  equation: 


Finally,  then,  division  of  (56)  by  (57)  gives 


Ai 

Ai 


-  <^)1/3 


-y 


1,0 


o.i 


(58) 


The  solutions  to  this  equation  are  the  eigen- 
energies  .  Note  that  (58)  replaces  the  simpler  condition 


Ai(x  ■  0)  *  0 


(59) 


used  in  space-quantization  theory  for  a  triangular  well 
bounded  by  an  infinite  wall  (?S  =  °°) .  In  our  problem  the 
solutions  to  Eq  (58)  differ  negligibly  from  those 
to  Eq.  (59) »  as  shown  by  Stern. 

The  tunneling  current  density,  corresponding  to 
the  pure  outgoing  wave  (52)  in  Region  II,  can  be  written  at 

-JL  1/3^  . 

Trm2  2 


J 


(6o) 
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The  coefficient  T  is  given  by  (56): 


T 


2 


A1  <-?!■<,> 

Bi2(-y2_0) 


2 

a 


(61) 


As  remarked  earlier,  the  wave  function  throughout 
Region  II  is  small  enough  that  we  can  neglect  it  in  normalizing 
the  total  wave  function.  Thus,  we  write  for  the  normalization: 


o 

2  2 

I  a  Ai  dx  =  1  ,  (62) 

j 

-00 

which  yields 


2 

a 


1 

Ai2(l+f) 


f 


Ai 

Ai' 


(63) 


where  the  argument  of  Ai  and  Ai '  is  -y  .  . 

1,0  p 

From  (58)  we  know  that  (Ai/Ai')  <  0.01.  Further, 

^  i .  — 

for  the  first  ten  energy  levels  E  ,  y  «.  -fo  •  Hence 

f  <  0.1,  and  can  be  neglected  in  (63). 

Using  the  various  results  above  in  (60)  we  finally 
obtain  for  the  tunneling  current: 


where 


J 


2 


mi  i  ^ 
“2  ^  ^ 


exp  -(FN) 


(64) 


(FN) 


U4  -  E) 
«2 


(65) 


and  the  no-ation  (FN)  has  been  used  because  the  exponent  ir 
(64)  is  Fowler-Nordheim-like  (see  Eqs.  (2)  and  (3)). 

It  is  convenient  to  re-write  (64)  in  the  form 


J 


4tt 

Z  3/2 


(V2 

“2 


E 

exp  -(FN) 


,(66) 
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where  Eq  is  the  lowest  space-quantized  energy-level  in  the  well: 


E 

o 


2/3 


( 2m^  j 


1/3 


Z 

o 


(67) 


ZQ  being  the  first  zero  of  the  Airy  function:  Ai(ZQ)  =  0,  and 
ZD  =  2.34  Here  we  have  used  (59)  as  a  suitable  replacement 
for  (58;. 

The  factor  EQ/h  in  the  pre-exponential  of  (66)  has 

the  dimension  of  frequency;  indeed  it  is  the  frequency  of  the 

periodic  motion  of  the  electron  in  the  well.  Physically  it 

corresponds  to  the  frequency  with  which  the  electron  hits 

the  potential  wall  6  in  Fig.  h2  and  has  accordingly  been 

called  "the  attempt-to-escape"  frequency  in  the  literature. 

The  pre-expcnentiai  factor  in  (66),  corresponding  to  tunneling 

from  a  bound  state  in  a  triangular  well,  replaces  the  pre- 

2 

exponential  factor  B£  of  Eqs .  (2)  and  (3),  corresponding  to 
tunneling  from  a  Sommerfeld  electron  gas,  in  which  case  one 
more  appropriately  considers  the  particle  current  incident 
on  the  barrier . 


C .  Application  of  the  Theory 

The  theory  of  Section  B  has  been  applied  to  the 
tunneling  experiments  on  the  Si-Si02  system  described  in 
Section  2.  These  experiments  were  all  done  on  (100)  Si,  for 
which  the  appropriate  (weighted)  electron  effective  mass  in  the 
space-quantization  problem  is  *  0.9.  The  electron  effective 
mass  in  the  oxide,  obtained  from  the  experiments  is  0.4  (see 
Eq.  (8)),  Eq  -  0.3  ev,  obtained  from  (67)  taking  ^  =  6  x  10^  V/cm 
(see  Fig  5).  The  corresponding  attempt-to-escape  frequency 
is  E^/'h  =  7  x  j.0"^/ sec .  Using  this  in  (66)  we  evaluate  the 
complete  pre-exponential  factor  as 


8  ** 
J  -  1.4  x  10  exp  -  (FN)  Amp/cm^ 


(68) 


The  corresponding  classical  Fowler-Noidheim  expression  (2)  and 


(3)  yields 


7  2 

J  -  l(  *  10  exp  -  (FN)  Amp/cm  .  (69) 

In  arriving  at  (68)  we  have  considered  tunneling 
only  from  the  ground  state  of  the  well.  Stern^”  has  estimated 
that  approximately  70#  oi  the  electrons  in  the  well  at  the 
interface  are  in  the  ground  state,  under  the  conditions  of 
our  experiment,  so  that  (68;  should  be  a  fair  estimate  for 
J,  even  allowing  that  the  remaining  30#  of  the  electrons  have 
a  higher  tunneling  probability  1  smaller  area  of  tunneling 
barrier).  The  small  image-iorce  corrections  to  (66)  are 
taken  up  in  a  more  complete  study  of  this  problem  to  be 
published  later. 

8.  Monte-Carlo  Studies  of  Hot-Electron  Behavior  in  Insulating 

Films  at  Hi^h  Fields  (S.  Baidyaroy  collaborating; 

A.  Background  Discussion 

In  the  First  Semiannual  Technical  Report  an  extensive 
discussion  of  our  Monte  Carlo  studies  of  hot  electrons  in  high 
fields  in  insulating  films  was  given.  Most  of  that  discussion 
was  concerned  with  constant -mean-f ree-path  problems,  with  both 
isotropic  and  anisotropic  scattering-  In  particular,  a  detailed 
study  was  made  of  one-dimensional  random-walk  simulations  of 
three-dimensional  problems .  A  beginning  was  made  on  the  study 
of  runaway  electrons  in  situations  where  the  mean  free  path 
is  energy-dependent  and  increases  with  energy.  In  this  report 
we  present  some  additional  information  on  one-dimensional 
simulations  of  anisotropic  scattering  problems  and  carry  the 
study  of  runaway  electrons  considerably  further. 

B.  Hew  Results 

In  the  First  Semiannual  Report  it  was  shown  that 
constant  mfp  (mean  free  path)  problems  with  anisotropy  could  be 
well  simulated  by  one-dimensional  random  walks  with  a  single 


anisotropy  parameter.  In  particular,  in  Fig.  28  of  that  report 

it  was  shown  that  the  computer  results  could  be  represented 

2 

by  the  simple  constant  mfp  scaling  law:  E  «  (FA)  /t  , 

o»v6 ,  ss  pn 

with  E  the  mean  energy  of  the  electron  emerging  from  the 

ave y  ss 

film  in  the  steady -state  (i.e.  independent  of  film  thickness), 

F  the  electric  field  intensity  in  the  film,  A  the  constant  mean 
free  path,  and  e  ^  the  optical  phonon  energy.  In  that  Fig. 
curves  were  plotted  for  f  *  0.2,  0.4,  0.5,  0.6  and  0.8,  f  being 
the  anisotropy  parameter  (ratio  of  forward-to-back  scattering 
probabilities).  A  mysterious  change  in  slope  appeared  in  the 
f  =  0.2  curve.  The  calculations  for  f  *  0.2  have  been  done  over, 
with  the  new  results  presented  in  Fig.  43.  Also  plotted  are 
new  results,  namely  the  curves  for  f  =  0.1  and  f  =  0.9  respec¬ 
tively,  and  the  old  results  for  f  =  0.4,  0.5,  0.6  and  0.8, 

Since  the  seven  plots  of  Fig.  43  are  all  linear,  the  relation¬ 
ship 

E  •  k  iHl-  ( 70) 

ave ,ss  f  e  .  * 

ph 

holds,  where  <  =  *(f)  is  an  f-dependent  constant  which  is 

determined  empirically,  namely  from  calculations  such  as 
exhibited  in  Fig.  43.  A  semi-long  plot  of  vb  f  is  g. ven 
in  Fig.  44,  with  points  available  in  the  range  of  0.1  _<  f 
0.99.  It  is  seen  that  in  the  'middle'  range,  0.1  _<  f  _<  0.8, 
a  simple  power-law  obtains: 


0.1  <  f  <  0,8  :  k  =  0.1  exp  4.22f  (71) 

Since  <  must  diverge  as  f  -*■  1  (pure  forward  scattering) 
this  relation  cannot  be  valid  for  f  close  to  1,  as  indeed  we 
note  in  Fig.  44. 

A  second  area  in  which  important  progress  has  been 
made  concerns  the  runaway  of  electrons  in  situations  where  the 
mean  free  path  is  energy-dependent  and  increases  as  the  energy 
increases.  Such  situations  were  examined  briefly  at  the  end 
of  the  First  Semiannual  Report,  in  Figs.  38- 4l  inclusive.  The 


_  no 


key  figure  was  Fig  38  which  is  reproduced  here,  with  several 

additicnaj.  curves  cn  it,  as  Fig.  45.  in  this  Fig.  the  average 

energy  E  m  the  emerging  electron  distribution  is  plotted  vs 

mfp  a.  Seen  on  the  Fig.  are  a  single  parabola  and  two  families 

01  straight  lines.  The  parabola  is  the  plot  of  E  8  0.87 

?  6 

i FA)  / e  ,  for  F  =  xO  V/cm.  and  £  8  0.1  ev ;  the  straight 

ph  pn 

lines  are  plots  cf  the  relations  a  =  \q  +  with 

r  =  0  i  ev  .  This  is  a  particularly  simple  form  of  energy- 
Ph 

dependence  of  the  mfp,  chosen  for  convenience.  The  parabola  is 

the  empirical  relationship  between  the  emerging,  steady-state 

average  energy  E  and  the  constant  ( energy— independent ) 

av  6  ,  s  s 

mfp  A,  assuming  isotropic  scattering,  previously  discussed  in 

the  First  Semiannual  Report.  The  two  families  of  straight 

o  ° 

lines  correspond  tc  a  *  25  A  and  a  *  40  A  respectively.  The 

o  o 

members  of  each  family  correspond  to  different  values  of  cu 
o  =  05,  0.8,  0-9,  1.0,  i.i,  1.2  and  2.0  for  the  'aq  =  25  A' 
family,  and  a  =  0  5  and  0 .8  for  the  'A  =  40  A'  family.  The 
nature  of  the  runaway  has  to  do  with  the  presence  or  absence 
of  intersections  of  the  straight  lines  with  the  parabola.  In 

O 


the  absence  of  intersections  (the  o  *  2.0  line  for  A  =  25  A 

o  u 

and  the  a  =  0.8  line  for  A  =  40  A)  there  is  strong  runaway, 

O  o 

as  seen  m  the  a  ®  1.0  curves,  for  A^  =  40  A,  in  Fig.  1+6,  and 
the  a  =  2.0  curves  in  Figs.  47  and  1+8  (as  well  as  in  Fig.  1+1 
of  the  First  Semiannual  Report).  For  the  case  that  there  are 
two  intersect  lens  of  the  straight  lines  with  the  parabolas 
the  situation  had  not  been  resolved  at  the  time  of  writing  of 
the  First  Semiannual  Report,  he  available  results  for  this 
situation  being  summarized  in  Fig.  39  of  that  Report.  A  more 
extensive  study  of  this  situation  has  led  to  a  considerable 
understanding  which  we  now  discuss. 

First  consider  the  two  a  =  0.5  curves  in  Fig.  1+6. 

From  Fig.  45  this  is  clearly  a  case  with  two  widely  separated 
intersections,  Yet  it  is  clearly  seen  in  Fig.  1+6  that  the  two 
a  =  0.5  curves  in  Fig.  46  do  not  converge  to  the  same  horizontal 
asymptote,  that,  in  fact,  the  curve  for  Eq  =  10.0  ev  actually 
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Fig.  45.  Graphical  Analysis  for  Runaway  Electrons 
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Fig.  U8.  Average  Energy  vs.  Film  Thickness  (A) 
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slowly  diverges  at  large  film  thicknesses.  The  reason  for  this 

interesting  behavior  is  that  for  a  batch  of  particles  entering 

at  the  energy  Ec  =  10.J  ev,  a  few  'escape'  from  the  well-behaved 

part  of  the  distribution.  ^These  few  'escaped'  electrons  are 

very  much  like  Shockley's  'lucky'  electrons).  A  particle, 

once  escaped  to  sufficiently  high  energy  (which  is  just  the 

second,  higher  i  '“rsection  point  in  Fig.  45)  is  almost  certain 

to  continue  gaining  energy  monotonically .  A  single  such 

particle,  in  a  finite  batch  of  particles,  is  enough  to  push 

the  average  energy  into  divergence  as  the  film  thickness  is 

taken  larger  and  larger.  The  curve  for  a  *  0.5  and  E^  =  0.01  ev 

has  not  a  single  escaped  electron  in  the  batch  of  200  particles, 

hence  the  convergence  of  a  wexl-defined  value  of  EfiVe  ss* 

The  one  or  two  escaped  particles  in  the  a  =  0.5,  Eq  =  10.0  ev 

tun  cause  the  slow  divergence.  Clearly  the  non-escaped  particles 

are  behaving  'propeny',  i.e  converging  as  expected  te.g., 

the  point  at  t  =  6000Aj.  Going  to  Fig.  47,  we  see  that  even 

for  E  *  0.01  ev,  lor  a  ac  low  as  0.8  there  are  escaped  particles, 
o 

even  though  there  are  two  well-separated  intersections  of 

o 

the  parabola  by  the  a  =  0.8,  aq  =  25A  straight  line  in  Fig.  45. 

For  X  =  25A  and  E  =10,0  ev  we  see  from  Fig.  48  that  there 
o  o 

are  escaped  particles  for  all  values  of  o  studied,  namely  down 
to  a  =  0.5.  Finally,  in  Fig.  49  we  compare  the  divergence  of 
the  a  *  0.9,  Xq  =  25A  curves  for  Eq  =  0.01  ev  and  Eq  =  4.0  ev. 

The  problem  of  'escaping'  electrons  is  examined 
more  closely  in  a  Technical  Report  which  will  be  issued  in  the 
coming  months.  In  particular,  the  actual  histograms  of 
emerging  particle  energies  wiil  be  presented  in  this  Report. 

These  histograms  confirm  in  detail  the  'escaping'  electron 
picture  implied  by  the  Eave  vs  t  (film  thickness)  curves 
discussed  above. 
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9.  Theoretical  Modeling  of  Local  Breakdown  IB .  Ridley  collaborating) 

A  model  for  electrical  breakdown  in  SiC>2  films  has 
been  constructed  along  the  following  lines,  illustrated  schematically: 


Electron  injection  via 

Motion  of  positive 

Fowler-Nordheim  tunneling 
from  projections  at  the 
interface 

Joule 

heating 

-i 

Activation  of 
positive  ions 

■3 

ions  to  cathode  and 
resultant  field 

enhancement 

±  '  = - 

Lateral  flow  of  heat 

'activates’  neighboring 
projections  into  similar 
cycles _ 


The  model  assumes  that  there  are  irregularities  at  the 

cathode  interface  which  produce  a  geometric  enhancement  of  the 

field  at  the  irregularity.  Fowler-Nordheim  tunneling  of 

electrons  into  the  oxide  at  the  irregularity  produces  a  locally 

high  current  density  which  produces  strong  Joulean  heating. 

Positive  ions  bound  in  the  oxide  volume  are  thermally  released 

and  drift,  in  very  high  field,  to  the  irregularity.  The  pile-up 

of  positive  ionic  charge  further  enhances  the  field  at  the 

irregularity,  yielding  a  positive  feedback  cycle  and  a  thermal 

runaway.  The  lateral  flow  of  heat  activates  neighboring 

projections  into  comparable  positive-feedback  cycles.  In  SiC>2 

the  postulated  ions  are  very  likely  Na  .  The  required  total 

sheet  densities  of  Na  are  on  the  order  of  10^"  to  10  ^/cm  and 

are  somewhat  dependent  on  details  of  the  model.  Such  densities 

appear  reasonable,  based  on  the  recent  work  of  Kuhn  and 

28 

Silversmith. 

The  details  of  the  above  model,  based  on  a  considerable 

amount  of  work,  will  be  presented  in  a  Technical  Report.  Here 

we  present  just  a  single  Figure,  namely  a  comparison  of  the 

prediction  of  the  theory  for  the  breakdown  field  dependence  on 

film  thickness, £  vs  d,  with  the  experimental  results  of  Osborn 
29  B 

and  Ormond  ,  Fig.  50. 
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Fig.  50.  Dependence  of  Breakdown  Field  on  Film  Thickness  (Al-SiOg-Ci  Structure) 
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